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Abstract 

Bushes of normal modes represent the exact excitations in the nonhnear physical 
systems with discrete symmetries [Physica D 117 (1998) 43]. The present paper 
is the continuation of our previous paper [Physica D 166 (2002) 208], where these 
dynamical objects of new type were discussed for the monoatomic nonlinear chains. 
Here, we develop a simple crystallographic method for finding bushes in nonlinear 
chains and investigate stability of one-dimensional and two-dimensional vibrational 
bushes for both FPU-a and FPU-/? models, in particular, of those revealed recently 
in [Physica D 175 (2003) 31]. 
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1 Introduction 



This paper expands upon our previous work [1] where an outline of the general 
theory of the bushes of modes [2,3,4] was presented in connection with the 
FPU-chain dynamics. Nevertheless, it seems reasonable to recapitulate here 
some basic notions and ideas. 
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1.1 The concept of bushes of modes 

Bushes of modes can be considered as a new type of exact excitations in 
nonlinear systems with discrete symmetries, such as molecules and crystals. 
The simplest way for introducing this concept can be described as follows. 

Let us have an A^-particle Hamiltonian system characterized by a symmetry 
group Go in its equilibrium state. Let us also assume that this system permits 
the harmonic approximation and, therefore, we can introduce a complete set 
of normal modes. Note that every normal mode possesses its own symmetry 
group Gj which is a subgroup of the group Gq. Let us now excite only one nor- 
mal mode at the initial instant by imposing the appropriate initial conditions. 
We call this mode by the term "root mode" . Normal modes are independent 
of each other in the harmonic approximation. However, if we take into account 
some anharmonic terms of the Hamiltonian, the excitation will transfer from 
the root mode to a number of other normal modes with zero amplitudes at the 
initial instant, so-called "secondary modes". Because of certain selection rules 
for the excitation transfer between modes of different symmetry, the number 
of the secondary modes can often be rather small. 

Definition. The complete collection of the root mode and all secondary modes, 
corresponding to it, forms a bush of normal modes. The number of modes in 
this collection is the dimension of the given bush. 

To avoid any misunderstanding at this point, let us note that one can consider 
normal modes, determined for a linear system, as a basis for a decomposition 
of different dynamical regimes in the nonlinear system (for more details, see 
below) . 

The energy of the initial excitation turns out to be trapped in the bush simply 
because of the above definition. The number of modes in the bush does not 
change in time while the amplitudes of these modes do change. 

For many systems, we can find one-dimensional, two-dimensional, three-dimensional 
and so on bushes of modes. Note that one-dimensional bushes can be treated 
as the similar nonlinear normal modes introduced by Rosenberg forty years 
ago [5]. 

Every bush possesses its own symmetry group G which is a subgroup of the 
symmetry group of the Hamiltonian of the considered system. In the above 
discussed case, when the bush is excited by perturbing the root mode, the 
symmetry group G of the bush coincides with that of this root mode. The 
symmetries of all other modes of the bush are greater than or equal to its 
full symmetry G. The development of the above ideas leads to the following 
important statement of the general theory [1,2,4]: 
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Proposition. Different nonlinear dynamical regimes of a physical system can be 
classified by subgroups of its symmetry group Gq (so-called "parent" group) ^ . 

It is important to emphasize that the considered bushes of modes are symmetry 
determined objects: the sets of their modes do not depend on the interparticle 
interactions in the physical system. Taking into account the specific character 
of these interactions can only reduce the dimension of the bush. 

1.2 Some properties of bushes of modes 

In general case, we must speak about bushes of symmetry- adapted modes 
rather than normal modes. Indeed, we consider the symmetry determined 
bushes whose finding does not need any information about the interactions 
in the system. The symmetry-adapted coordinates are the basis vectors of 
the irreducible representations (irreps) of the parent symmetry group Go and 
they can be found by group-theoretical methods only, without any informa- 
tion about the interactions. In contrast, the normal coordinates, used for con- 
struction of normal modes, must be obtained, generally, by diagonalization of 
the Hamiltonian matrix and, therefore, the interparticle interactions must be 
taken into account. 

In the geometrical sense, a bush of modes is an invariant manifold singled 
out by its symmetry and decomposed into the basis vectors of the irreducible 
representations of the appropriate parent symmetry group. 

This decomposition plays a very important role in our approach. Indeed, dif- 
ferent irreps describe the transformation properties of the dynamical variables 
corresponding to different physical characteristics of the considered system. 
For example, some irreps are active in the infrared or Raman experiments, 
while others are nonactive in optics, but play the essential role in neutrono- 
graphics, etc. 

In this paper, we discuss bushes of vibrational modes describing the time- 
dependent displacements of the particles of nonlinear chains from their equi- 
librium positions. Speaking about the symmetry group of a single mode or of 
the bush as integral object, we mean the specific symmetry of the patterns of 
the instantaneous displacements of all particles which correspond to the mode 
or to the bush. Considering the symmetry group of the displacement pattern 
and the collection of the normal modes contained in a given bush, we deal 
with the geometrical aspect of the bush. This aspect is independent of the 
interparticle interactions in the chain. 

^ We prefer to use this term because the considered theory is valid not only for the 
Hamiltonian systems, but also for the dissipative systems. 
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On the other hand, we deal with the dynamical aspect of the bush when 
consider the differential equations describing the time-dependence of the am- 
plitudes of the bush modes. In this sense, the given bush represents a reduced 
dynamical system whose dimension is, in general, less (often, considerably 
less!) than the dimension of the original physical system. Obviously, the dy- 
namical aspect of the bush, in contrast to the geometrical aspect, depends 
essentially on the interparticle interactions. 

The modes of the bush are coupled by the so-called "force interactions" , while 
their coupling to all other modes is brought about by "parametric interac- 
tions" [4]. The latter interactions can lead to a loss of stability of the bush, 
if the amplitudes of its modes become sufficiently large [4,1,6], because of the 
phenomenon similar to the well-known parametric resonance. In this case, the 
spontaneous breaking of the symmetry of the system's vibrational state takes 
place and, as a consequence, the given bush transforms into a bush of larger 
dimension. Precisely this mechanism of stability loss by the bushes of the FPU 
chains is discussed in this paper. The other mechanisms of the loss of bush 
stability will be considered elsewhere. 

1.3 Some history 

The general concept of bushes of modes for physical systems with discrete 
symmetry was introduced in [2]. Actually, this concept appeared as a gener- 
alization of the notion of the complete condensate of primary and secondary 
order parameters, whose theory was developed in our papers devoted to phase 
transitions in crystals [7,8]. The discussion of properties of bushes of modes, 
the group-theoretical methods for their construction, some theorems about 
the bush structure can be found in [2,3,4]. The bushes of modes for different 
classes of mechanical systems with point and space symmetry were investi- 
gated in [3,4,6,9,10,11]. In particular, let us mention the investigation of 

- all possible "irreducible bushes" and symmetry determined similar nonlinear 
normal modes for the mechanical systems with any of the 230 space groups [9]; 

- low-dimensional bushes for all point groups of crystallographic symme- 
try [10]; 

- bushes of vibrational modes for the fuUerene Cqq ( "buckyball" structure) [11] ; 

- bushes of modes and their stability for the octahedral structure with Lennard- 
Jones potential [6] . 

A review of the group-theoretical methods and the appropriate computer algo- 
rithms for treating the condensates of order parameters and bushes of modes 
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can be found in [12]. A remarkable computer program ISOTROPY by Stokes and 
Hatch, capable to deal with many group-theoretical methods of crystal physics 
and, in particular, with bushes of modes, is now available on the Internet as 
free software [13]. 

In [14], Poggi and Ruffo revealed some exact solutions for the FPU-5 chain 
dynamics, by analyzing the appropriate dynamical equations, and investigated 
their stability. These solutions correspond to certain subsets of normal modes 
(so-called "subsets of I type") "where energy remains trapped for suitable 
initial conditions" . 

In [1], we showed that these subsets are the simplest cases of the bushes 
of modes for the monoatomic nonlinear chains and that they can be found 
without taking into account the information about the specific interactions 
in the FPU-/3 model. We discussed there the symmetry-determined bushes of 
modes with respect to the translational group T, in detail, and the bushes 
with respect to the dihedral group D, partially. Also, in the paper [1], the 
stability of the former bushes for the FPU-o; chain was studied. 

In a recent paper [15], Bob Rink developed a method for finding all possi- 
ble symmetric invariant manifolds for the monoatomic chains. In spite of the 
fact that this paper was carried out independently of our approach, based 
on the concept of bushes of modes, and that it is written in a more rigor- 
ous mathematical style, the main idea and the group-theoretical method are 
very similar to those in our previous papers. Actually, Bob Rink found all 
symmetry-determined bushes of modes for the monoatomic chains with re- 
spect to the dihedral symmetry group. Moreover, he revealed some bushes 
for the FPU-/? chain which are brought about by the additional symmetry 
of this mechanical system connected with the evenness of its potential (here- 
after we call them "additional bushes"). However, let us note that the general 
idea of the classification of the symmetric manifolds (bushes of modes) by 
the subgroups of the symmetry group of the Hamiltonian was not emphasized 
explicitly in this paper. 

Finally, we would like to comment on two papers by Shinohara [16], devoted 
to finding "type I subsets of modes" (bushes, in our terminology) in nonlinear 
chains with fixed cndpoints. The author docs not use any symmetry-related 
method and prefers to analyze the specific structure of the dynamical equa- 
tions (practically, in the manner similar to that of Poggi and Ruffo in [14]). 
Moreover, he writes that the "approach, which refies on the direct analysis of 
the mode-coupling coefficients, is much simpler" than the one based on the 
group-theoretical methods. In contrast to this opinion, we think that the sit- 
uation is quite the opposite! In particular, we have tried to demonstrate in 
the present paper that the symmetry-related methods are much simpler and 
much more transparent than those based on investigating the structure of the 



5 



dynamical equations. This simplicity becomes especially apparent when one 
deals with physical objects more complex than monoatomic chains, such as 
molecules and crystals characterized by nontrivial symmetry groups. 

Note that comparing our results and those from [16], one must take into ac- 
count that we use periodic boundary conditions, while Shinohara discusses 
the nonlinear chains with fixed endpoints (the former boundary conditions 
are more general, since the latter conditions can be considered as their special 
case) . We would also like to note that the stabihty problem for the bushes of 
modes (invariant manifolds or "modes subsets of I type") is not considered 
in [16]. On the other hand, this problem seems to be important for treating 
the induction phenomenon discussed in those papers. 

1.4 The structure of the present paper 

In Section 2, we consider a simple crystallographic method for obtaining 
bushes of modes in the space of atomic displacements which differ from that 
in our previous paper [1] and in the paper [15] by Bob Rink. A major advan- 
tage of this method is its remarkable clarity. Actually, we obtain the invariant 
manifolds corresponding to the subgroups of the dihedral group and then 
decompose them into the normal modes of the monoatomic chain. Here we 
restrict ourselves to the discussion of one- and two-dimensional bushes only 
(about other bushes see [17]). 

The dynamical description of the bushes of normal modes is presented briefly 
in Section 3. 

The general approach to analyzing the parametric stability of the bushes of 
normal modes is discussed in Section 4. 

Sec. 5 is devoted to investigation of the stabihty of the one-dimensional (ID) 
bushes in the FPU-a and FPU-/3 chains. Only one ID bush, B[a^,i], can be 
obtained for the monoatomic chain with an even number of particles ^ , if the 
translational group T is considered as the appropriate parent group of this 
mechanical system. The stability of the bush B[a^, i] was studied for the FPU- 
a and FPU-/3 chains in the papers [1] and [14], respectively.^ However, we 
obtain some additional one-dimensional bushes considering the dihedral group 
D as the parent symmetry group [1,17]. Moreover, as was already mentioned, 
there exist still more ID bushes for the FPU-/3 model because of the evenness 
of its potential [15]. Therefore, it is interesting to compare the regions of 

^ It was denoted by the symbol B[2a] in the paper [1]. 

^ The stability of this bush, consisting of only one so-called vr-mode, was also dis- 
cussed in [18]. 
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stability of all these bushes in both FPU-a and FPU-/? models. We perform 
such a comparison in this section. An exceptional case is represented by the 
bush B[a^, ai] for the FPU-o; chain for which the threshold of stability loss is 
exactly equal to zero (this result can be proved analytically). 

We discuss here the stability of ID bushes not only for the FPU chains with 
the finite number of particles {N), but for the continuum limit {N — > oo), as 
well. 

In Section 6 we discuss the shape of the stability regions of the two-dimensional 
(2D) bushes for the FPU-a and FPU-/9 chains. Indeed, thresholds of the loss 
of stability of the ID bushes depend only on their energy, while the similar 
thresholds for multi-dimensional bushes depend on the complete set of the 
initial conditions for the dynamical equations of the bush. We demonstrate, 
with the aid of the appropriate plots, how nontrivial the stability regions of 
the 2D bushes for the FPU chains can be. The dependence of these regions on 
the total number of particles in the chains is partially discussed. 

In Conclusion (Sec. 7), the obtained results are summarized and some new 
directions for the bush theory are outlined. 



2 Simple crystallographic method for finding bushes of vibrational 
modes in the monoatomic chains 

We consider the longitudinal vibrations of A^-particle monoatomic chains with 
periodic boundary conditions. Let the A^-dimensional vector 

X{t)^{xi{t),X2{t),...,XN{t)} 

describe all the displacements, Xi{t), of the individual particles (atoms) from 
their equilibrium states at the instant t. In accordance with the above-mentioned 
periodic boundary conditions, we demand for arbitrary time t: 

XN+l{t) = Xi{t). (1) 

In the equilibrium state, a given chain is invariant under the action of the oper- 
ator a which shifts the chain by the lattice spacing a. This operator generates 
the translational group 

r = {^,a,a^...,a^-^}, = ^, (2) 

where E is the identity element and N is the order of the cyclic group T. The 
operator d induces the cyclic permutation of all particles of the chain and. 
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therefore, it acts on the "configuration vector" X{t) as follows: 

dX{t) = d{xi{t),X2{t), . . . ,XN-l{t),XN{t)} = {xN{t),Xi{t),X2{t), . . .,XN-l{t)}. 

The full symmetry group of the monoatomic chain contains also the inversion 
i, with respect to the center of the chain, which acts on the vector X{t) in 
the following manner: 

iX{t) = i{xi{t),X2{t), XN-l{t),XN{t)} = {-XN{t), -XN-l{t), -X2{t), -Xi{t)}. 

The complete set of all products d'^i of the pure translations a*^ (/c = 1, 2, . . . , N— 
1) with the inversion i forms the so-called dihedral group D which can be writ- 
ten as a direct sum of two cosets T and T ■ i: 



The dihedral group is non-Abelian group induced by two generators (a and i) 
with the following generating relations 



Let us consider the geometrical interpretation of the symmetry elements of 
the dihedral group D. For simplification, we will use the refiection a in the 
mirror plane orthogonal to the chain instead of the inversion i. Indeed, the 
actions of a and i on the vectors of the one-dimensional space are equivalent 
to each other {ax — ix — —x). Then we can use the well-known theorem of 
crystallography: "Any mirror plane and the translation A orthogonal to it 
generate a new mirror plane, parallel to the former plane and displaced away 
from it by A/2". 

Let the number of particles in the chain (N) be an even integer. Then it is 
easy to check that the products d'^a = d^i {k = ±1, ±2, ±3, . . . ) generate 
the family of parallel planes depicted by vertical lines"' in Fig. 1. For odd k, 
these planes pass through the atoms, while for the even k they pass between 
atoms, exactly at the middle of the neighboring particles. Note that, for odd 
N, everything will be vice versa: for even k, the planes pass through the atoms 
and for the add k - between them. In Fig. 1, we show a part of the family 
of the above mentioned planes near the center of the chain. It can be verified 
that the symmetry elements d'^a = d^i {k = 1,2,3,...) are located to the 
right of the center of the chain, while the elements ad'' = id'' are located to 
the left from it ^ . 

^ For simplicity, the hats above the operators are dropped in all our pictures. 



(3) 



d^ — E, i'^ — E, id — d ^i. 



(4) 



^ Note that a'^i = id '' = id 



.N-k 
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Fig. 1. Symmetry elements of the dihedral group D near the center of the 
monoatomic chain. 

Let us now consider the subgroups Gj of the dihedral group D, keeping in 
mind that every subgroup Gj C D induces the certain bush of modes. 

Each group Gj contains its own translational subgroup Tj C T", where T is the 
above discussed full translational group (2). If is divisible by 4 (for example, 
we consider below the case N = 12) there exists the subgroup T4 = [a^] of the 
group T = [a]. Note that in square brackets we write down only generators of 
the considered group, while the complete set of group elements is written in 
curly brackets (see, for example, Eqs. (2)). 

If a vibrational state of the chain possesses the symmetry group T4 = [a^] = 
{E, a^, a^, . . . , d^~'^}, displacements of the atoms, which are apart by the dis- 
tance 4a from each other in the equilibrium state, turn out to be equal. There- 
fore, for the case N — 12, we can write the following atomic displacement 
pattern: 

X {t) = {Xi {t) , X2 {t) , X3 (t) , X4 {t) I Xi {t) , X2 {t) , {t) , X4 {t) I Xi {t) , X2 {t) ,Xs{t),X4 (t) } , 

(5) 

where Xi{t) {i — 1,2,3,4) are arbitrary functions of time. The pattern (5) 
determines the bush corresponding to the symmetry group T4 = [a^]. Thus, 
the complete set of the atomic displacements can be divided into N/4 (in our 
case, N/4: = 3) identical subsets, and each of them determines the so-called 
"Extended Primitive Cell" (EPC) ^ . For the bush (5), the EPC contains four 
atoms, and the vibrational state of the whole chain is described by three such 
EPC. In other words, the EPC for the vibrational state with the symmetry 
group T4 = [a^] (its size is equal to 4a) is four times larger than the primitive 
cell (a) of the chain in the equilibrium state. 

It is essential that some symmetry elements of the dihedral group D disappear 
as a result of the symmetry reduction D = [a, i] ^ T4 = [a'^]. In our case, there 
are no nontrivial symmetry elements of the group D inside the EPC and, as a 



^ This term is used in the theory of phase transitions in crystals. 
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Fig. 2. Vibrational bush corresponding to the subgroup G = [a'^, i]. 

consequence, there are no restrictions on the displacements belonging to one 
and the same extended cell (this fact is obvious from Eq. (5)). 

There are four other subgroups of the dihedral group D to which the same 



translational subgroup T4 



correspond: 



[d [d ,di], [d ,d i], [d ,d i]. 



(6) 



These subgroups possess two generators, namely, the common translational 
generator and different inversion elements d'^i {k — 0, 1, 2, 3). These inver- 
sion elements differ from each other by the position of the center of inversion 
(see Fig. 1 where they are depicted by the vertical lines representing the mirror 
planes) . 

Note that subgroups [a^, d'^i] with A: > 3 are equivalent to those from the 
list (6), because the second generator a^i can be multiplied by d^'^, represent- 
ing the inverse element with respect to the first generator d'^. Thus, there exist 
only five subgroups of the dihedral group (with N mod 4 = 0) constructed 
on the basis of the translational group T4 = [a^], namely, this group and the 
four groups from the hst (6). 

Now, let us consider the bushes corresponding to the subgroups (6). 
1. The subgroup [a^,i] consists of the following six elements: 

E,d'^,i,d'^i,d^i ^ id'^. (7) 



The diagram of this subgroup, similar to that of the full dihedral group D in 
Fig. 1, is depicted in Fig. 2. 

Two symmetry elements (ia^ and a^i), situated at the boundaries of the EPC, 
do not produce any additional restrictions on the atomic displacements inside 
this EPC. Indeed, they connect with each other the atomic displacements in 
the adjacent extended primitive cells, while the equivalence of the displace- 
ments of the atoms with the same positions in the adjacent EPC is guaranteed 
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by the translation a which is the first generator of the considered subgroup 
[a^^]. 

On the other hand, the inversion i at the center of the EPC (Fig. 2), surviving 
in the subgroup [a^, i], demands the displacements of the atoms, symmetrically 
situated with respect to this center, to be equal in value, but opposite in sign: 

Xi{t) = -X^t), X2{t) = -X3{t). (8) 

These relations between the atomic displacements are shown in Fig. 2 by the 
appropriate arrows. As a consequence, the atomic displacements, correspond- 
ing to the bush B[a'',i] with the symmetry group [a^,i], can be written (for 
N — 12) in the following form: 

X{t) = {xi{t),X2{t), -X2{t), -Xi{t) I Xi{t),X2{t), -X2{t), -Xi{t) \ Xi{t),X2{t), -X2{t), -Xi{t)} 

(9) 

Since Eqs. (8) hold for an arbitrary moment t, we will drop the argument t 
and write down the vibrational bush in the X -space (configuration space) as 
follows: 

B[a^,i] = |a;i,a;2, -a;2, -xi 
Thus, we point out the atomic displacements in only one EPC. 
2. For the subgroup 

[a^, di] = {E, a^, d^, di, d^i = id'^ , d^i = id^}, (11) 

we obtain the diagram depicted in Fig. 3. In this diagram, one can see two 
symmetry elements, id^ and di, which pass, respectively, through the first and 
third atoms of the represented EPC. Therefore, the displacements of these 
atoms must satisfy the relations xi{t) = —xi{t), 2:3 (t) = —X3{t). In turn, 
this means that Xi{t) = 0, X3{t) = 0, i.e. the atoms located on the inversion 
elements (which are equivalent, in the one-dimensional space, to the mirror 
planes!) must be immovable. 

On the other hand, the symmetry element di is situated in the middle be- 
tween the second and fourth atoms of the EPC and, therefore, their possible 
displacements must satisfy the equation X2{t) = —X4{t). Thus, we have the 
following displacement pattern for the bush 3(0"^, di] 

X{t) = {0, x{t), 0, -x{t) I 0, x{t), 0, -x{t) I 0, x{t), 0, -x{t)} (12) 



(10) 
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Fig. 3. Vibrational bush corresponding to the subgroup G = [a^, ai]. 
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Fig. 4. Vibrational bush corresponding to the subgroup G = [a^, a^i]. 

with only one independent dynamical variable x(t). We can write this one- 
dimensional vibrational bush in the abbreviated form: 

B[a^, ai] — |0, x, 0, —x\. 

3. The diagram of the symmetry elements of the subgroup 

[a^ a^i] = {E, a^ a^ a^i, a^^i = ia\ a}^i = ia"} 

is depicted in Fig. 4. There are two elements, id'^ and d'^i surviving in the 
transition from the full dihedral group D = [d,i] to the considered subgroup 
[a^,a^i]. These elements are located between two first and two last atoms 
of the depicted EPC. Therefore, xi{t) — —X2{t), xs{t) — —X4{t) and the 
corresponding displacement pattern is 

X (t) = {Xi {t),-Xi {t) ,X2{t),-X2{t)\Xi {t) , -Xi {t) , X2 {t) , -X2 {t) \ Xi {t) , -Xi (t) , X2 {t) , -Xg (t) } . 

(13) 

There are two degrees of freedom in (13), Xi(t) and X2{t). Thus, the two- 
dimensional bush with the symmetry group [a^, a^i] can be written as 

B[a^,a^i] = |a;i, -a;i,a;2, -X2I. 
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Fig. 5. Vibrational bush corresponding to the subgroup G = [a'^,d^i]. 
4. The diagram of the symmetry elements for the last subgroup of the hst (6) 

[d'^, d^i] = {E, a^, a®, d^i, d^i = ia^, d^^i = id] 

is depicted in Fig. 5. We see from this diagram that X2{t) = 0, X4{t) = 0, 
xi{t) — —X3{t), and, therefore, 

X{t) = {x{t), 0, -x{t), I x{t), 0, -x{t), I x{t), 0, -a;(t), 0}. (14) 

The corresponding one- dimensional bush can be written as follows 

B[a'^, d'^i] = \x, 0, —x, 0|. 

We have just considered all the subgroups of the dihedral group D, which 
contain the common translational group T4 = [a^], and the corresponding vi- 
brational bushes. In the same way, we can find the subgroups of the group 
D corresponding to the arbitrary translational group = [a"*] ^ and then 
obtain the bushes induced by these subgroups [17]. We will not discuss this 
method any further in the present paper, because all the bushes of modes (sym- 
metric invariant manifolds) for the dihedral group were recently obtained by 
Bob Rink in [15]. Nevertheless, let us note that the above described crystal- 
lographic method for finding bushes of vibrational modes seems to be more 
straightforward than that used in [15]. Actually, the former method repre- 
sents a variant of our general method of the splitting of the space group orbits 
(SO) [12,8], adapted for the sufficiently simple case of the dihedral group D. 
The SO method was developed for finding the so-called complete condensate 
of primary and secondary order parameters in the framework of the theory of 
phase transitions in crystals. It was applied to such nontrivial space groups as 
01 [8], 01 [19], etc. 

The notion of the complete condensate of order parameters is equivalent to 
the notion of the bush of modes in the geometrical (only geometrical!) sense. 
Three different methods for finding the condensate of order parameters were 

^ Such subgroup exists in the case N mod m = 0. 
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considered in [7,12,8]: the "direct method", which is similar to that by Bob 
Rink [15], the above mentioned SO method and the method based on the 
concept of the irreducible representations of the appropriate symmetry groups. 
The last method was used, in particular, for treating bushes of vibrational 
modes in nonlinear chains in the paper [1] . 

Now let us consider the notion of the bush equivalence. There are two pairs of 
equivalent bushes of modes induced by the subgroups from the list (6). Indeed, 
shifting the displacement pattern (12) of the bush B[a^, di] by the lattice 
spacing a, we obtain the displacement pattern (14) of the bush B[a^,a'^i]^. 
The pattern (9) of the bush B[a^,i] can be transformed into the pattern (13) 
of the bush B[d'^, a^i] by the same shifting. 

The importance of the notion of bush equivalence is brought about by the 
fact that equivalent bushes represent idenUcal dynamical systems: they are 
described by equivalent equations of motions (see below). It can be proved, in 
the general case, that the conjugate subgroups of the parent symmetry group 
induce equivalent bushes ^ . Remember that two subgroups Gi and G2 of the 
same group G (Gi C G, 6*2 C G) are called conjugate to each other (Gi ~ G2), 
if there exists at least one element g^ of the group G which transfers Gi into 
G2 by the transformation 

G2^g^'Grg^ (goeG). (15) 

It can be easy checked the following conjugations: 







"-4 -2^' 




"^4 






d^, i 




a ,a t 


5 


a , at 




a ,a t 



For example, using Eqs. (4), we obtain the relation 

[a^, di] — d~^ ■ [a^, d^i] ■ a, 

which proves the equivalence of the subgroups [a*^,ai] and [a'^,a"^i], if we take 
into account the definition (15) with go = d. Thus, the both above considered 
one-dimensional bushes B[a'^,ai] and B[a*^,a^i] turn out to be equivalent to 
each other, as well as the both two-dimensional bushes, B[a^, i] and B[a^, a^ij. 
Moreover, it can be proved that for any even m and A; = 0, 1, 2, . . . , m — 1, all 
bushes of the form B[a™, a^i] with k of the same evenness are equivalent. For 
odd m, all the bushes of this form are equivalent. 

Let us summarize the above said about the construction of the vibrational 
bushes for nonlinear monoatomic chains. 

^ We can imagine that our chain represents a ring, because of the periodic boundary 
conditions. 

^ Sometimes, we call equivalent bushes by the term "dynamical domains". The 
term "domain" is borrowed from the theory of phase transition in crystals. 
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- Each subgroup of the dihedral group singles out a certain bush of vibrational 
modes. 



- Every subgroup is determined by one or two generators and contains the 
translational subgroup — [a''] (it is trivial, To, for the subgroups [a/^i], 
m = 0,l,2,...,Ar-l). 

- The translational subgroup = [a''] determines the extended primitive cell 
(EPC) for the vibrational state of the chain. 

- The symmetry elements a™i, belonging to the second coset of the dihedral 
group (see Eq. (3)) and surviving at the symmetry lowering D = [d,i] 
[a'^, a"^i], produce certain restrictions on the displacements inside the EPC. 

- The conjugate subgroups induce the equivalent bushes of modes. The num- 
ber of the equivalent bushes of a given type is equal to the index of the 
corresponding subgroup in the dihedral group. 

In addition, note that there exist subgroups to which no vibrational bushes 
correspond. In our case, the subgroup [a^,ai] does not induce the vibrational 
bush, because the EPC, corresponding to this subgroup, contains only two 
atoms through which the inversion elements, ia and di, pass and, therefore, 
these atoms cannot move. 

The bushes obtained with the aid of the described method represent invariant 
manifolds determined in the configuration space. Let us now consider the 
bushes in the modal space. 

We consider the normal modes for the monoatomic chain in the form used 
in [14]: 



■ . f27ik \ f27ik ^ 

sm n + cos n 

\ N \ N , 



n = l,2,...,7V . (17) 



Here the subscript k refers to the mode, while the subscript n refers to the 

atom. The modes ■j/'^ (A; = 0, 1, 2, . . . , — 1) form an orthonormal basis in the 
modal space and we can decompose the set of the atomic displacements X{t), 
corresponding to a given bush, in this basis 

N-l 

Xit) = J2 Mi)'^k (18) 

k=0 

For example, we obtain the following decompositions for the bushes B[a'^,i] 
and B[a^, a^i]: 



The index of the subgroup G of the group Gq is equal to the ratio ||Go||/||G| 
where ||Go|| and ||G|| are the orders of these two groups, Gq and G. 
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X{B[a^,i]) = {xi{t),X2{t), -Xi{t), -X2{t)\xi{t),X2{t), -Xi{t), -X2{t)\xi{t),X2{t), -xi{t), -X2{t)} 



= ll{t)ll)N/2 + J^{t)ll}3N/4, (19) 

X(B[a^ aH])^il{t)^N/2 + m'^N/A- (20) 

Only vectors i/'Af/2) ''/'JV/4 ^"^^ from the complete basis (17) contribute 

to these two-dimensional bushes: 

11,^/2 = ^(-1, 1, -1, 1, -1, 1, -1, 1, -1, 1, -1, 1), (21) 
^^^^74 = ^(1, -1, -1, 1, 1, -1, -1, 1, 1, -1, -1, 1), (22) 
^3iv/4 = ^(-1, -1, 1, 1, -1, -1, 1, 1, -1, -1, 1, 1). (23) 



Remember that these bushes are equivalent to each other, i.e. they represent 
the so-called "dynamical domains". For the bush B[a^,i] we can find the fol- 
lowing relations between the old dynamical variables Xi(t), X2{t) (relating to 
the configuration space) and the new dynamical variables //(i), fit) (relating 
to the modal space): 

Mt) = -^[a:i(t)-X2(t)], 

v{t) = -^[x,{t)+X2{t)]. 

Thus, each of the above bushes consists of two modes. One of these modes is 
the root mode {iI^^n/a for the bush B[a^,i] and iI^n/a for the bush B[a^,a?i\), 
while the other mode '4>n/2 is the secondary mode. Indeed, we have already 
discussed (see Introduction) that the symmetry of the secondary modes must 
be higher or equal to the symmetry of the root mode. In our case, as it can 
be seen from Eqs. (22), the translational symmetry of the mode iI^n/a is 
(acting by this element on (22) we obtain the same displacement pattern), 
while the translational symmetry of the mode ■0iv/2 is a^, which is two times 
higher than that of iI)n/a (see (21)-(23)) . 

One-dimensional (ID) and two-dimensional (2D) vibrational bushes for the 
FPU chains in configuration and modal spaces can be seen in Tables 1,2. Note 
that in these tables, as well as in Tables 3,4,5 we list one-dimensional bushes 
completely, but two-dimensional bushes partially: only those 2D bushes are 
given whose stability properties are studied in the present paper (see Sec. 6). 
The complete list of 2D bushes in the FPU chains can be found in Appendix. 



Note that the full symmetry of the modes ■07v/4 ^-^id ipN/2 ^^e [a"^, a^i] and [a?, i], 
respectively. 



(24) 
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Table 1 

Representation of ID and 2D bushes of vibrational modes for the FPU chains in 
the configuration space (classification by the dihedral group D) 



Bush 


l~)iRnlfi,rpTnPTit DRttprn 


Dimension 




|x, — x\ 


1 


B[a^] 


Ti To To * 


2 


B[a^,i] 


'r* f) — T 


1 


Bffl^ ni] 

-l—f \ IA/ ■ 1X6 


f) n" — 


1 




|x, — X, 0| 


1 


B[a^,z] 


\x^ Xo — Xo — X^ 1 


2 


Bfa^ a^i] 




2 


B[a^, ai] 


0, X, 0, —X 


1 


B[a4,a3i] 


X, 0, —X, 0| 


1 


B[a^a^] 


|0,Xl,X2,0,-X2,-Xl| 


2 




|X1,0, — Xl,X2,0, — X2I 


2 




.r2. 0. -.r2. -.I'l. 


2 



* Here X3 is not an independent variable: its value must be determined by the 
condition of immobility of the mass center of the chain, xi + X2 + X3 = 0. 

Table 2 

Representation of ID and 2D bushes of vibrational modes for the FPU chains in the 
modal space (classification by the dihedral group D). Here £1 = sin (^) = ^Hi^/^zH^ 

V2(v/3+l) _ _ %/2 
4 ) =-3 — 2 



£2 = cos 



(n) 



Bush 


Representation in modal space 


B[a2,i] 




B[a^,i] 

B[a^,ai] 

B[a^,a'^i] 


fi^li^N/S + £2lp2N/3) 
^{£2tpN/3 + £lV'2Ar/3) 
l'£3(V'iV/3 - V'2iV/3) 


B[a^^] 
B[a\a^ 


Z^lV'3Ar/4 - ^2tl^N/2 
Z^l'0iV/4 + ^'2'0Ar/2 


B[a'^, ai] 
B[a^,a3i] 


'^e3('0Af/4 +'03JV/4) 
'^e3(''/'Ar/4 -'03JV/4) 


B[a^ai] 

B[a6,o3i] 

B[a^a5^] 


Z^l(eit/>Ar/6 + £2lp5N/6) + 1^2 ( - ^2 -^Af/S " £l^2Ar/3) 
iyi{£2lpN/6 + eiV'SiV/e) + '^2(£l'0JV/3 + £2ll^2N/3) 
Z^l£3(V'Ar/6 - V'SiV/e) + t'2£3(V'Ar/3 - ■02Ar/3) 
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Table 3 



Additional ID and 2D bushes, for the FPU-/? chain (these bushes were obtained 
in [15]) 



Bush 


Dim. 




T?pnrPRPTitatimi in tbp moHal snapp 


B[a^,iu] 




x, —2a;, a; 


vieotbNi-^ — ei'02/V/3) 


B[a^, aiu] 


1 










T 7! — 2t 




B[a^, ml 

L ' J 


1 


1 X ^ 3C ^ X ^ X 1 




13 [f7 ^ . f/^ /'('/] 




.r. — — .r 


' -iA/ 1 


B[a^, ai, a^u] 




|0, X, X, 0, — x, —x\ 




B[a^,a^i,a^u\ 


1 


1 x J ^ x J X ^ ^ X 1 


^{£2lpN/6 + ei^5Af/6) 






X, X, 0, —X, —X, 


J^£3(.ll^N/6 - ^57V/6) 


B[a*^, aiu] 


2 


Xl , Xo, Xs, Xol** 


V^e-^ilbN/A — '^■^N/a) + Z^2'0Af/O 






|X2,X1,X2,X3|** 




B[a\a^u] 


2 


\xi,X2, — Xl, — X2I 




B[a^,m] 




Ixi Xo Xq X^j Xo Xi I** 

1^1? ^0? ^05 "^il 


7/1 f Pn lh\Tla + £1 l/jc; Ar /^! ) + fo fSo lb\Tlr> — £^ ibn J\.T lr> ) 
■'iV'^zV^iV/D ^ i'roly /oj ' ly^ 2.^^ N j 6 y'r IN 1 i) 


B[a^,a2m] 


2 


|xi,Xi,X2,X3,X3,X2|** 


Z^l£3(V'7V/6 - ■05Ar/6) + J^2£3{tJ^N/3 + tl^2N/3) 


B[a^, a'^iu] 




|xi,X2,X2,Xi,X3,X3|** 


UlieilpN/Q + e2'05Ar/6) + '^2(ei'«/'A'/3 " £2ll^2N/3) 


B[a^, i, a^u] 




|xi,X2,Xi, -Xl, -X2, -Xil 


Uii-eilpN/Q + e2^5Af/6) + '^2'0Ar/2 


B[a^,a'^i,a^u\ 


2 


|X1, -Xl,X2, -Xl,Xl, -X2I 


Z^iesl'j/'TV/e + V'5Af/6) - '^2^Ar/2 


B[a^, a'^i, a^u] 




|xi,X2,-X2,-Xi,-X2,X2| 


Z^l(£2'«/'iV/6 - £lV'5JV/6) + i^2'0JV/2 



** See footnote to Table 1. 



3 Dynamical description of the bushes of modes 

Up to this point, we've been considering bushes as geometrical objects. In 
the present section, we address their dynamical aspect. For this purpose, one 
would need to obtain the equations of motion for the appropriate dynamical 
variables (in the configuration space or in the modal space). It can be done 
with the aid of the Lagrange method. This method is very useful for the me- 
chanical systems with some constraints on the natural dynamical variables 
(for example, on the displacements of the individual particles). Note that 
we deal now with precisely this case, and the additional constraints on the 
dynamical variables brought about by the symmetry causes. Indeed, the dis- 
placement pattern X(t), corresponding to a given bush B[G], is determined by 
the condition of its invariance under the action of the bush symmetry group G: 
GX{t) — X{t). Namely this condition permits us to introduce the generahzed 
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coordinates which represent dynamical variables of the reduced Hamiltonian 
system associated with the considered bush of modes. 

We have already used the Lagrange method for the derivation of the bush 
dynamical equations in the paper [6] devoted to nonlinear vibrations of the 
octahedral mechanical systems with Lennard- Jones potential. Now we apply 
this method to the FPU chains. 

The Hamiltonian of these nonlinear chains can be written as follows: 

N . N N 

H^T + V^-J2xI + -Y: {Xn+1 - Xnf + ^ E - Xnf ■ (25) 

^ n=l ^ n=l P n=l 

Here p = 3, 7 = « for the FPU-a model, and p = 4, 7 = for the FPU-/9 
model . T and V are the kinetic and potential energy, respectively. We also 
assume the periodic boundary conditions (1) to be valid. Let us consider the 
set of atomic displacements corresponding to the two-dimensional bush B[a^, i] 



X{t) = {x, y, -y, -x\x, y, -y, -x\x, y, -y, -x\--- }, (26) 

where we rename Xi{t) and X2{t) from Eq. (9) as x{t) and y{t), respectively. 
Substituting the atomic displacements from Eq. (26) into the Hamiltonian (25) 
corresponding to the FPU-a chain we obtain the following expressions for the 
kinetic energy T and the potential energy V: 

T^^{x' + f), (27) 



V = j{3x^ - 2xy + 3y^) + ^{x^ + x^y - V - y')- (28) 

These equations are valid for the arbitrary FPU-o; chain for which mod 4 = 
0. The size of the extended primitive cell for the vibrational state (26) is equal 
to 4a and, therefore, we can restrict ourselves calculating the energies T and 
V, by summing over only one EPC. 

The Lagrange equations can be written as 

0, (29) 




dxj 



where L = T — V . The subscript j is used here for numbering the dynamical 
variables (in our case, xi{t) = x{t), X2{t) = y{t)). Using Eqs. (27) and (28), we 



The constants a and j3 can be removed from the Hamiltonian (25) by a trivial 
scaling, but we keep them for the clarity of consideration. 
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obtain from (29) the following dynamical equations for the considered bush 
B[a^^]: 



X + (3x - y) + (xi'ix^ + Ixy - y^) = 0, 

V yj \ y y J , ^^^^ 

y+{3y-x) + a{x'^ - 2xy - 3y'^) = 0. 

Let us emphasize that these equations do not depend on the number N of the 
particles in the chain (but the relation N mod 4 = must hold) . 

Eqs. (30) are written in terms of the atomic displacements x{t) and y{t). From 
these equations, it is easy to obtain the dynamical equations for the bush in 
terms of the normal modes /i(t) and i^it). Using the relations (24) between the 
old and new variables, we find the following equations for the bush B[a^, i] in 
the modal space 



/i + 4/x--=i/2 = 0, (31) 



u + 2iy =aiy = 0. (32) 



The Hamiltonian for the bush B[a*^,?], considered as a two-dimensional dy- 
namical system, can be written in the modal space as follows: 

H[a\ i] = lifi' + ^') + + u') - (33) 



We already emphasized that identical dynamical equations correspond to 
equivalent bushes of modes. Moreover, the dynamical equations of the given 
type are often associated with many bushes of very different symmetry groups 
(these equations differ from each other by numerical values of their pertinent 
parameters only). Because of this reason, we can say that such bushes be- 
long to the same class of "dynamical universality" [2,3]. As an example, let 
us point out that all one-dimensional bushes for the FPU-/? models belong to 
one and the same class of dynamical universality. Indeed, the dynamics of all 
these bushes is described by the Duffing equation jl -\- Aji -\- B ji^ — 0, but with 
different values of the constants A and B. 

Dynamical equations for all one-dimensional and two-dimensional bushes of 
vibrational modes for the FPU-chains arc presented in Tables 4,5. In the first 
column, we present the complete sets of equivalent bushes, to which the same 
dynamical equations (see column 2 and 3) correspond. 

Table 5 describes the "additional" bushes , which were found in [15]. They 
They are additional with respect to the bushes of normal modes associated with 
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Table 4 



Dynamical equations for ID and 2D bushes in the FPU chains, induced by the 
subgroups of the dihedral group 



Bush 


Equations for the FPU-o; chain 


Equations for the FPU-/3 chain 


B[a'^,i] 


+ = 


iV 


B[a^,i] 
B[a^, ai] 
B[a^,a^i] 


^ VJV 


z iV 


B [a^ , ai] 


i> + 2z/ = 


/; + 2z/ = -4^=1/3 


B[a^] 


V2 + 3z/2 = ^^(-i^? + 2l/iZ/2 + I'o) 


^i + 3z/i = -f +^^2) 


B[a^,i] 
B[a\a^i] 


i>2+ 41/2 = -4^1/2 


i>l + 2z/i = -8;^Z/i(z/2 + 3j/2) 

i>2 + 4z/2 = -8|:i/2(3z/2 + 2i/2) 


B[a^,ai] 

B[a^,a^i] 

B[a^,a^i] 


i^i + z^i = -\/6^z/iz/2 

^2 + 3^2 = ^^(-^^1^ + 31/2) 


^1 + 1^1 = -i|''^iK^ + 3i/|) 

^2 + 3z/2 = -|^I^2(z^f + 3l/|) 



are brought about by the evenness of the potential of the FPU-/5 model and 
the additional symmetry operator u is used in the description of these bushes. 
This operator changes the signs of all atomic displacements: uX{t) = —X{t). 
Note that the operator u is contained in the description of the bushes only in 
combinations with some other symmetry operators, i, or their products. 



4 Stability of bushes of normal modes 

The term "stability" is often used in various senses. Let us explain what we 
mean by this term in the present paper. 

In general case, stability of the bushes of modes was discussed in [2,4], while 
their stability in the FPU-o; chain is considered in [1]. As well as in these 
papers, we will discuss here the stability of a given bush of normal modes 
with respect of its interactions with the modes which do not belong to this 

the dihedral group. 

There is an essential difference between the interactions of the modes which belong 
and which do not belong to a given bush: we speak about "force interaction" in the 
former case and about 'parametric interaction" in the last case [4] . 
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Table 5 

Dynamical equations for the additional ID and 2D bushes, for the FPU-/3 chain 



Bush 


Equations 


B[a'^,iu] 
B[a^, aiu] 




B[a^, m] 
B[a^,a^m] 


i; + 2i/ = -8|^z^3 


B[a^, ai, a'^n] 
B[a^, a-^i, a^u] 
B[a^,a^i,a^u] 




B[a^, am] 
B[a^, a^iu] 


i>i + 2i.i = -44z/iK + 6z/|) 
i;2 + 4l/2 = -^^vi{Zv\ + 2f|) 


B[a^a2^x] 


i>i + 21^1 = -8^z/f 
£>2 + 4i^2 = -S^z^l 


B[a^^^t] 
B[a^,a^m] 
B[a^, a'^iu] 


+ = -fl^z^iK^ + 9f|) 


B[a^, i, a^u] 
B[a^,a^i,a^u] 


,72 + 47/2 = -|^(\/2/4^ + \2vlv2 + 167yf) 



bush. Let us illustrate this idea with the following example. 

For the FPU-a chain, the two-dimensional bush B[a^, i] is described by Eqs. (31), (32). 
These equations allow a solution of the form 

/i(t)7^0, z/(^)=0. (34) 

(The vibrational regime of this type can be excited by imposing the appro- 
priate initial conditions: /x(to) = /^o 7^ 0, /i(to) = 0, vito) = 0, z>(to) = 0). 
Substitution of the solution (34) into (31) produces the dynamical equation 
of the one-dimensional bush B[a^, i] (see Table 5) consisting of only one mode 
/x(t) 

/i + 4/i = (35) 



22 



with the trivial solution 



IJ,{t) = fiQ cos(2i). 



(36) 



Substituting (36) into Eq. (32), we obtain 



i>+ 2- 



cos(2i) i/ = 0. 



(37) 



This equation can be easily transformed into the standard form of the Mathieu 
equation 



On the other hand, there exist domains of stable and unstable motion of the 
Mathieu equation (38) in the a — q plane of its intrinsic parameters, a and q. 
The one-dimensional bush B[a^, i] is stable for the sufficiently small amplitudes 
/io of the mode /u(t), but it becomes unstable as the result of the increasing 
of this amplitude. This phenomenon, similar to the well-known parametric 
resonance, takes place for those /xq which get into the domains of unstable 
motion of the Mathieu-type equation (37). The loss of stabihty of the dynam- 
ical regime (34) (the bush B[a^,i]) manifests itself in the appearance of the 
mode v{t) which was identically equal to zero for the vibrational state (34). 
As a result, the dimension of the original one-dimensional bush B[a^,i] in- 
creases and this bush transforms into the two-dimensional bush B[a'^,i]. This 
transformation is accompanied by the breaking of the symmetry of the vi- 
brational state (the symmetry of the bush B[a^,i] is twice larger than that 
of the bush B[a^,i]). In general case, we consider a given bush as a stable 
dynamical object, if the complete collection of its modes (and, therefore, its 
dimension) does not change in time. All other modes in the system, according 
to the definition of the bush as the full collection of active modes, possess zero 
amplitudes, they are "sleeping" modes. If we increase intensity of bush vibra- 
tions, some sleeping modes, because of the parametric interactions with the 
active modes, can lose their stability and become excited. In this situation, we 
speak about the loss of stability of the original bush, since the dimension of 
the vibrational state (the number of active modes) becomes larger, while its 
symmetry becomes lower. As a consequence of the stability loss, the original 
bush transforms into another bush of higher dimension. 

As it was described above, the loss of stability of the bush B[a^, i] with respect 
to its transformation into the bush B[a^,i] can be investigated with the aid 
of the Mathieu equation. But if our nonlinear chain contain particles, we 
must analyze the stability of the bush B[a^, i] not only relatively to the mode 

We can consider the initial phase in the solution (36) to be equal to zero [1]. 



i>-F [a - 2qcos{2t)\v = 0. 



(38) 
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of the bush B[a^,i], but with respect to all the other modes, too. This 
can be done by the Floquet method. 

The stabihty of the bush B[a^,i] ^® in the FPU-a chain for arbitrary even 
N was investigated by us in [1]. Let us recall some points of this investiga- 
tion which are needed for analyzing the stability of the other one-dimensional 
bushes. 

In the modal space, the dynamical equations of the FPU-a chain can be 
written as follows (see Eqs. (27) from the paper [1]): 



, (39) 



••2 r , \ /r. • 27r . 27r .\ 

uj + UjUj = --j= ^ Uk ^{uj+k + i^-j-k) (2 sm —k + sm —j j 

+ (i/j_fc - i/_j+fc) \2 sm —k - sm — j 

u;| = 4sin2(^), j = 0,l,2,---,7V-l. 

Here Vjify are the time-dependent coefficients in front of the basis vectors 
t/^j (17) in the decomposition (18). For brevity, we call Vj{t^ by the term 
"mode" , in spite of the fact that this term, rigorously speaking, corresponds 
to the product v^{t)t\}y The numbers of all modes Vj in (39) are assumed to be 
reduced to the interval l<i<A^ — Iby adding ±.N ^ since these numbers are 
defined modulo N . Hereafter, the mode vq is excluded from the consideration 
of the vibrational bushes because it corresponds to the motion of the chain as 
a whole. 

The bush B[a^,i] consists of only one mode VN/^it) and looks as 

B[a^i] = i/iv/2(t)^/'jv/2, (40) 

where ■j/'jv/2 is determined by Eq. (21). The dynamical equation for the mode 
^N/2{t) can be obtained from (39) supposing that all the other modes are equal 
to zero: ^^{t) = 0, j 7^ y. This equation reads 

+ l^N/2^N/2 = 0. (41) 

Its solution can be written in the form 



PN/2{t) = Acos{2t) (42) 

(the initial phase in this solution, as well as in Eq. (36), may be chosen equal to 
zero [1]). Linearizing the system (39) near the exact solution (40) with respect 

In the paper [1] , this bush was denoted by the symbol B [2a] . 
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of all the modes uj (j y), we obtain the following approximate equations 
.. 2 Sa . /2nj\ 

uj + Uj ^ --^sm\^—j vn/2T^n/2-3 , (43) 



9 9/7r?\ N N 

^J = 4sin2(^^j, j = l,2,.-.,--l,- + l,-.- ,7V-l. 

The system (43) splits into a number of subsystems containing one or two 
equations only. Indeed, the mode Vj in (43) is connected with the mode i/j, 
where j — y — J , and vice versa, the mode Vj is connected with = ^j- 

Therefore, wc have for j = 1, 2, . . . , y — 1 (see Eqs. (43) the following pairs of 
linearized equations which are independent from all the other equations: 



i)j +4sm [§jiyj = -7 sin (^j z/at/s-j cos(2r), ^^^^ 
^N/2-j +4cos2 (21 j i'N/2-j = -7 sin (^) i^j cos(2r), 

SaA N 



Here we replaced i'N/2{t) with 74cos(2t) from Eq. (42). For j — . . . ,N — 

1, we obtain equations equivalent to Eqs. (44) [1]. 

If j = ^ — j and, therefore, j — Eqs. (44) reduce to a pair of identical 
equations of the form 

i>jv/4 + 2z/iv/4 = -Wn/4 cos(2r), (45) 

which is the Mathieu equation. Actually, this is the above considered case, 
describing the transformation of the bush B[a^, i] with the mode fi{t) = i'N/2{t) 
into the two-dimensional bush B[a^, i] containing the modes //(i) = i'N/2{t) and 

It is convenient to rewrite Eqs. (44) as follows: 

X + A:sm^{k)x — 7sin(2/i;)y cos(2t) ^^^^ 
y + 4cos^(A;)y = 7sin(2A;)xcos(2T), 

where k = x{t) = Vjit) and y{t) = i'N/2-jif)- Thus, studying the loss of 
stabihty of the bush B[a^, i], brought about by its interactions with the modes 
(and i^zn/a)-, is reduced to analyzing the Mathieu equation (45), while the 
loss of its stability with respect to all other modes reduces to analyzing the 
equations (46). 

Eqs. (46) are linear differential equations with periodic coefficients and, there- 
fore, they can be studied with the aid of the Floquet theory. The system (46) 
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is very interesting, and wc wonld like to cite the following fragment of the text 
of our paper [1], devoted to its properties: 

"Our computation of the multiplicators for the system (46) as eigenvalues 
of the monodromic matrix reveals a surprising fact! Indeed, it seems that 
the critical value 7c of the constant 7 from (44), corresponding to the loss 
of stability of the bush B[a^,i], must depend on the mode number j, since 
the coefficients of these equations depend explicitly on j^. However, this is 
not true. We found that 7c does not depend on j^, at least up to 10~^, and 
coincides with that of the Mathieu equation (45) : 

7c = 2.42332. (47) 

Then from Eq. (44) we obtain 

aA^ = 0.30292yiV. (48) 

This nontrivial fact means that the original bush B[a^, i] loses its stability with 
respect to all modes Uj (j 7^ y) simultaneously, i.e. for the same value of 
the amplitude of the mode i'N/2{t)- In other words, all modes of the A^-particle 
FPU-a chain are excited parametrically because of interactions with B[a^,i] 
when aA reaches its critical value aAc — 0.30292-\/]V, and, therefore, the bush 
B[a^,i] transforms at once into the bush B[a^] of trivial symmetry". 

One may think that the above-mentioned fact means that there exists such 
transformation of Eqs. (44) which removes the dependence of these equations 
on the parameter (or the dependence on k of Eqs. (46)). Nevertheless, this 
doesn't turn out to be correct. Indeed, in Fig. 7, wc show (in white color) 
not only the first stability domain for the different modes for < < vr 
{k = ^), but also the second stability domain. The former domain looks as 
a white horizontal band near the A;-axis, while the latter domain looks as two 
white triangles above this band (the vertical axis corresponds to the value A 
of the amplitude of the i'N/2{t) mode of the bush B[a^,i]). This means that 
there is no dependence of the excitation conditions for the modes Uj (j 7^ y) 
on the number j for the first stability domain, but such a dependence is very 
obvious for the second stability domain. 

Since we use similar stability diagrams for all one- dimensional bushes, let us 
describe their structure in more detail. In Fig. 7, each point (A, k) determines 
a certain value of the root mode amplitude A of the bush B[a^, i] and a certain 
value k — j^, which is connected with the index j of a fixed mode of the FPU-a 
chain. The black points {A, k) correspond to the case where the mode j = k— 
becomes excited because of its parametric interaction with the mode of the 
bush B[a^,i]. The white color denotes the opposite case: the corresponding 
mode j, being zero at the initial instant, continues to be zero in spite of its 
interaction with the considered bush. Dashed vertical lines correspond to the 
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modes (j = 1, 2, . . . , 12) for the FPU chain with = 12 particles. Considering 
stabiUty diagrams, we must make one remark concerning normaUzation of the 
normal modes. In Eq. (17), there is the normalizing factor which leads 

to the relation IV'fcl^ = 1- Because of this factor, the absolute values of the 
individual atomic displacements corresponding to a given mode t/'fe tend to 
zero for N ^ oo. Such normalization of the normal modes is not convenient 
for studying the bush stability. Indeed, the critical vahic of atomic displace- 
ments, associated with the loss of stability, say, for the bush B[a^,i] is equal 
to Xc = 0.30292/q; (this result can be easy obtained from Eq. (48) (see [1]). 
This value Xc does not depend on the number N of particles in the FPU 
chain. Because of this reason, depicting stability diagrams of one-dimensional 
and two-dimensional bushes in Figs. 28-33, we use different normalization of 
the normal modes than that in Eq. (17), namely, we simply reject the factor 
which presents in Eq. (17). Then, the maximal value of the atomic dis- 
placements corresponding to every normal mode t/j^ turns out to be equal to 
V2. 



We have just discussed the stability of the bush B[a^,i]. Now we intend to 
consider the stabihty of all the other one-dimensional bushes in the FPU-a 
and FPU-/3 chains. There are three such bushes for every monoatomic non- 
linear chain with dihedral symmetry group, but for the FPU-/3 model, there 
exist three additional one-dimensional bushes because of the evenness of its 
potential [15]. For investigating stability of these bushes (three for the FPU-a 
chain and six for the FPU-/3 chain), we shall also use the Floquet method. 



There are two difficulties in implementing of this method for one-dimensional 
bushes in the general compared to that for the above considered bush 

B[a^,i]. Indeed, the dynamical equations for the bushes B[a^,i] and B[a*^,ai] 
in the FPU-a model are reduced to the harmonic oscillator equation (with 
different frequencies), while the dynamical equations for all the other bushes 
turn out to be nonlinear. As a consequence, the period T of the oscillations, 
described by the bush of such a type, depends on the amplitude Hq of these 
oscillations (T = T(/io)). Therefore, we must construct the monodromy matrix 
by integration of the linearized dynamical equations, for the considered chain, 
over the period T{fj,o). This period must be found for every value /Iq before 
obtaining the monodromy matrix. 



The second difficulty is that the linearized dynamical system (near the ap- 
propriate exact solution to the given bush) can be split into the subsystems 
whose dimension (3 or 4) can be higher than that for the case of the bush 
B[a^, i]. As a result, we must deal, respectively, with 6 x 6 or 8 x 8 monodromy 
matrices. 



27 



5 Stability of one-dimensional bushes for the FPU chains 



We will consider the stability of the bushes for the FPU-o; and FPU-/3 models 
separately 



5.1 FPU-a chain 



The exact dynamical equations for the FPU-a chain in the modal space are 
given by Eqs. (39). 

We want to linearize these equations in the vicinity of the solution described 
by a given one-dimensional bush. There are three one-dimensional bushes as- 
sociated with the subgroups of the dihedral group of the monoatomic chains 
with the appropriate divisibility properties of the number of the particles 
(see Table 2): B[a^, i] for N mod 2 = 0, B[a^ i] for N mod 3 = 0, B[a^ at] for 
N mod 4 = 0. The following root modes correspond, respectively, to them: 
i'{t)ij)N/2, J^(t)[eiil^N/3 + £2'4^2N/3], i^{t)b/^N/4: + il^3N/4]- Notc that the root modes 
for the bushes B[a'^, i] and B[a^, di] are the certain superpositions of the conju- 
gate modes with time-independent coefficients Si — sin (^^^ , £2 = cos ) 



Let us consider the bush B[a^, a^i] . Linearizing the FPU-« dynamical equa- 
tions near the solutions Xq^I) = i'{t)ss[ipN/s — 'tp2N/3] of this bush £3 = ^ , 
we must consider the amplitudes i'N/3{t) = £3i^{t) and z^2JV/3(^) = —£3^{t) in 
Eqs. (39) to be much greater than those of all other modes or, more precisely, 
iyj{t) with J 7^ ^ must be treated as the infinitesimal values. Then we re- 
tain in Eqs. (39) only linear terms with respect to the modes ^'j(t) (j ^ y, ^) 
which do not contribute to the given bush. The result of this computational 
procedure can be written as follows 
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v{t) 



{l + ^sin 
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(49) 



J 



1,2, 



N -I. 



According to [1], we call two real modes and if^N-j by the term "conjugate" 
since normal modes and (j>N-j, corresponding to them in the complex form, are 
complex conjugate. The frequencies of the conjugate normal modes are equal. 
^® This bush, being equivalent to the bush B[a^, i], turns out to be more suitable for 
the next calculations. 
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It can be easy seen from this form of linearized dynamical system that, for 
every fixed j, the modes Uj, I'N/s+j and z^2JV/3+j, being coupled by Eq. (49), are 
decoupled from all other modes. Therefore, we obtain for j = 1, 2, ■ ■ ■ , y — 1 
the following 3x3 systems of ordinary differential equations with periodic 
coefficients proportional to the function ^{t) : 



u{t) [(l - 2 cos + f )) u.^, + (l - 2 sin + |)) u 

+ 4+/f = -^^(i) [(l - 2 cos + f )) z., + (l + 2 cos 2^ 
+ '^k+/^+,- = -^^(t) [(l - 2 sin + f )) i^j + (l + 2 cos 



3 

I Z/2JV 



N J ^^+3_ 



(50) 

In turn, the function v{t) is the solution to the dynamical equation 



3v^ 



a 



i) + 3u — ■;=U 



,2 



of the considered bush B[a^,a^i]. Thus, the stability of the bush B[a^,a^i] 
can be investigated by the Floquet method applied to the systems (50) for 
j = l,2,--- ,f -1. 

Similarly to the above said, we can obtain the 4x4 systems of differential 
equations for investigation of stability of the bush B[a^,ai]. All systems of 
such a type, which must subjected to the further Floquet analysis, are given 
for the FPU-a chain in Table 8. 



5.2 FPU- (3 chain 



We have obtained the linearized dynamical equations for one-dimensional 
bushes in the FPU-a chain by linearizing the exact nonlinear equations (39) 
defined in the modal space. This method, namely, obtaining the exact equa- 
tions in the modal space with their subsequent linearization, is not the sim- 
plest method for our purpose. Indeed, we can, firstly, linearize the dynamical 
equations in the configuration space and only then transform the obtained 
equations to the modal space. In this manner, we can obtain the linearized 
equations in the modal space for different nonlinear chains, in particular, for 
the FPU-/9 chain. Moreover, we want to find the dynamical equations lin- 
earized near the solution 

X^{t) = y{t)c, (51) 

which describes the exact nonlinear regime corresponding to the considered 
one-dimensional bush. Here c is a constant vector which determines the di- 
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rections of displacements of all particles in the chain, while uit) is a time- 
dependent function satisfied the dynamical equation of a given ID bush 

i) + uPv ^ F{v). (52) 
Here F{y) is a certain nonlinear function. 
Let us then suppose that 

X{t) = Xo(t) + 5{t) (53) 

with 5{t) — {5i(t),-- - ,5Ar(i)} being an infinitesimal vector, and substitute 
X{t) in this form into the dynamical equations of the nonlinear chain in the 
configuration space. For the chain with dynamical equations 

Xk = (xfc+i+a:fe_i-2a;fc)+7 [{xk+i - Xk)'^ - {xk - Xk-i)""] , A; = 1, 2, • • • , A^, 

(54) 

and with periodic boundary conditions x^^iit) = Xi{t), XQ{t) = xis[{t), after 
linearization with respect to S{t), we obtain 

6 = B{t)d. (55) 

Here time-periodic matrix B{t) is almost three-diagonal (see Fig. 6) and can 
be determined by the equations 

= -2 - 7(i) [(c,+i - c^r-' + {cj - cj-^r-'] , 

^'.J-i = 1 + 7^(9 -c,-i)—i, 
6,,,+i = l + 7(t)(c,+i-c,)— 1, 

where ^{t) — ^mu'^~^{t). 

The vector Eq. (55) represents the hnearized dynamical equations of the con- 
sidered chain (54) in the configuration space. Now, we must transform these 
equations to the modal space. Let 

N-l 

Ht) = E ^.WV'.' (56) 

where the basis vectors are determined by Eq. (17). Substituting d{t) 
from Eq. (56) into Eq. (55) and taking into account that these basis vectors 
are orthonormal, one can obtain the final dynamical equations of the chain, 
linearized near the solution (51) to a given one-dimensional bush. It must be 
noted, that this transformation, being simple in principle, is practically very 
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For the FPU-a and FPU-/3 chains, m is equal to 2 and 3, respectively. 
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Fig. 6. The structure of the matrix B from Eq. 55 (for N = 8). 

cumbersome. We performed it with the aid of the mathematical packet MAPLE. 
The results of this procedure are hsted for all one-dimensional bushes in the 
FPU-a and FPU-/3 chains in Table 8 and Table 9, respectively. 

As can be seen from Table 9, the full systems of linearized differential equations 
corresponding to different one-dimensional bushes in the FPU-/3 chain split 
into certain subsystems whose dimensionalities are equal to 1, 2, 3.^*^ 

These systems contain time-periodic coefficients proportional to i^'^{t), where 
is the solution to the Duffing equation 

iy + uP'v — . 

The coefficients cu^ and 7 are, naturally, different for different bushes. 

Using Tables 8 and 9, we can investigate the stability of all one- dimensional 
bushes, with the aid of the Floquet method, for the FPU-o; and FPU-/? models, 
respectively. 

In our bush stability diagrams (see Figs. 7-15), we depict regions of sta- 
bility (white color) and instability (black color) for individual modes j = 
1, 2, 3, . . . , — 1 of the FPU chains. The loss of stability by any mode, which 
does not belong to a given bush, leads to the loss of stability by the entire 
bush. The fact that we can study stability of individual modes allows us to 
investigate the bush stability not only for the chains with finite number of 
particles, but also for the case of thermodynamical (continuum) limit — > 00 
(see a discussion in the next section). 

"^^ Let us note that splitting of the general dynamical system, linearized near a 
given bush, into independent subsystems of small dimensionalities, is brought about 
by symmetry-related causes. Some general results concerning such splitting will be 
considered elsewhere. 
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Practically, we study the eigenvalues of the monodromy matrices for a certain 
subsets of sleeping modes. We choose as many points on the horizontal fc-axis 
(A; = as it is necessary for obtaining sufficiently accurate picture. In some 
cases, the density of the analyzed A;-points corresponds to ~ 10^. Usually, 
we recognize instability of a given set of sleeping modes, if there are eigenvalues 
of the corresponding monodromy matrix with modules exceeding 1 by 10~^. 

In Table 10, for the case = 12, we present the thresholds of the loss of stabil- 
ity for all one-dimensional bushes in the FPU-a and FPU-/3 chains. In the sta- 
bility diagrams (Figs. 7-15), all possible modes for this case (j = 1, 2, 3, ... , 12) 
are depicted by vertical dashed lines. Using these diagrams, one can under- 
stand, in some sense, the cause of the given threshold values of bush stability 
for the FPU chains with = 12 (see below). 

5.3 Diagrams of stability for one- dimensional bushes in the FPU-a chain 

First of all, we would like to emphasize, that studying the stability of bushes 
of normal modes, we can restrict ourselves to only one specimen from the set 
of "dynamical domains" which represent equivalent bushes (see the previous 
section). 

Now we consider the stabihty diagrams for the bushes in the FPU-a chain. 

5.3.1 BushB[a^,i] 

The stability diagram for this bush (see Fig. 7) has already been discussed. 
Let us recall that in all our diagrams, for one-dimensional bushes, the stability 
regions are white, while unstable regions are black. 

5.3.2 BushB[a^,i] 

In the configuration and modal spaces this bush looks, respectively, as follows: 

X{t) = {x{t), 0, -x{t)\x{t), 0, -x{t)\ ■■■} = iy{t) {eiViv/s + £2V'27v/3} • (57) 

The stability diagram for the bush B[a'^,i] is presented in Fig. 8. The domain 
depicted in grey color corresponds to the disruption of the FPU-a chain 
(atoms escape from the potential wells in which the vibrational motion can 
occur). Thus, we must distinguish between the stability of the chain itself and 

Obviously, this phenomenon is impossible for the FPU-/3 chain. 
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the stability of a given vibrational regime, i.e., of a bush of modes, in this 
chain. 



Let us imagine a horizontal line in Fig. 8 corresponding to a given value A of 
the bush root mode. This line can partially pass through the white and par- 
tially through the black regions of the considered stability diagram. Actually, 
the parts of this line belonging to the black (unstable) region represent the 
set of modes which are excited because of the interactions with the bush root 
mode at the fixed level of its amplitude {A). If A is very small, our horizontal 
line A = const intersects the narrow unstable regions near the mode numbers 
i = 0,f,f^ and N. 

Let us suppose that A'^ = 12 and A « 1. Then it is easy to see from Fig. 8 
that only two modes, corresponding to j — ^ 3 — excited in 

our mechanical system This fact is consistent with the stability of the 
bush B[a'^, i] for small values of its root mode amplitude. Indeed, according to 
Eq. (57), this bush contains only two (conjugate) modes i^N/sit) and i'2N/3{t) 
with the certain relation between their amplitudes: 

J^N/3{t) = £liy{t), 1^2N/3{t) = £2J^{t) (58) 

Moreover, we can find the threshold of the stability loss of the bush B[a^,i] 
directly from the diagram in Fig. 8, taking into account that all the possible 
modes for the case = 12 are marked by vertical dashed lines. This threshold 
is the minimal vertical distance from the horizontal coordinate axis to the 
black (unstable) region for j = 1, 2, 3, 5, 6, 7, 9, 10, 11. The stability thresholds 
for all one- dimensional bushes in the FPU-a and FPU-/3 chains for A^ = 12 
are given in Table 10. From this table, we find that the stability threshold for 
the bush B[a^,i] (the critical value Xc of the atomic displacements) is equal 
to Xc = 0.203. The corresponding critical values of the root mode amplitude 
{Ac) and the energy {Ec) are equal to Ac = 0.166, Ec = 0.047. 

Let us now increase the value A of the root mode amplitude of the bush 
B[a^,i]. Such an increase leads to the broadening of the unstable regions near 
the j = 0, y, which intersects the line A = const in Fig. 8. This fact 

does not affect the stability of the bush B[a'^,-i] for A^ = 12, as long as there 
are no modes uj with j = 1, 2, • • • , 11 (except the modes un/s and i^2iv/3) which 
fall into the above-mentioned broadening unstable intervals. 

On the other hand, if A^ ^ oo (continuum limit), the mode density, i.e. the 
density of vertical dashed lines in Fig. 8, increases and there appear narrow 
groups of modes near the values j — 0, y,^,A^ {k = 0,|,^,7r). This case 
corresponds to the onset of the modulation instability of the considered bush of 



Note that the mode vq (and the equivalent to it mode 1^12) must be excluded 
because of the momentum conservation law [1]. 
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modes. Note that we then have a bush of narrow wave packets which graduahy 
collapses on a long time scale Thus, we conclude that the bush B[a"^,i], 
being stable for small number of particles in the chain (iV), for a fixed value of 
its root mode amplitude, becomes unstable in the continuum limit — > oo. 



5.3.3 BushB[a^,ai] 

The one-dimensional bush B[a*^, di] can be written as follows 

a/2 r 

X{t) = {0,x{t),0,-x{t)\0,x{t),0,-x{t)\- ■ ■} = ^Vit) + V'STVA 

(59) 

Thus, it represents a superposition of two conjugate modes "^^n/a and iI^zn/a 
with equal amplitudes 

^N/iit) = l'2,N/i{t) = ^^^(^)- (60) 



According to the Table 10, the threshold of the stabihty loss for the bush 
B[a^, di] in the FPU-a chain, in contrast to the similar thresholds for all other 
bushes, is equal to zero. In particular, the loss stability threshold (LST) for 
the same bush in the FPU-/3 chain is equal to 1.161 (for the case N — 12). 
Let us discuss this unusual case in more detail. 

In Fig. 9, we see that besides the modes VN/i{t) and v^N/iit) contributing to 
B[a^,ai], the mode i'N/2(t) is also excited for arbitrarily small values of the 
root mode amplitude of the considered bush. On the other hand, all the results 
presented in Table 10 and in Figs. 7-15 are obtained numerically by the Flo- 
quet method and, partially, by a straightforward integration of the dynamical 
equations of the FPU-chains. Since numerical results are only approximate, 
it is desirable to prove analytically that the LST for the bush B[a^, di] in the 
FPU-a chain is indeed equal to zero exactly. 

In the general case, the presence of three modes I'N/iit), i'N/2{t), i^m/iit)) in 
the dynamical regime, means that we have the three-dimensional bush B[a^] 
(see [17]) with 

X{t) = l'N/A{t)ll)N/A + l'N/2{t)'tl)N/2 + i^3Af/4 (t)'03Ar/4- (61) 

Therefore, one can suppose that the initially excited one-dimensional bush 
We will discuss this phenomenon elsewhere. 

This is the upper boundary of the first stabihty domain of the considered bush. 
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B[a^, ai] transforms spontaneously into the three-dimensional bush B[a^]: 



B[a*,ai] ^B[a^]. (62) 

This transformation is accompanied by the lowering of the symmetry, be- 
cause of the loss of the invertion element di. In the dynamical picture of the 
above transformation, the new mode i'N/2{t) appears and the specific relation 
between the modes z^Ar/4(t) and J^37v/4(^), presented in the bush B[a^,ai], is 
destroyed: 



Thus, we must now prove that the bush B[a^, a,i] loses its stability with respect 
to the interaction exactly with the mode i'n/2- So, let us exclude all the other 
modes of the FPU-a chain from our consideration and focus only on the 
transformation (62). Obviously, if we prove that the stability threshold of the 
bush B[a^,ai], for the transformation (62), is zero, then it is also zero with 
respect to the interactions with all modes of the chain. 

Because of the above said, let us consider the dynamical equations of the 
three-dimensional bush B[a^] which can be written as follows 

ill + 2/11 = -^fiifj,2, 

A2 + 4/.2 = -f (/x?-/.i), (63) 



1X3 + 2/^3 = ^//2/^3- 



N 
N' 



Here we renamed the modes ^^/^{t), i^N/2{t) and u^N/iit) as /ii(t), H2{t) and 
Ii3{t), respectively. 

In the dynamical regime corresponding to the bush B[a^,ai] only one degree 
of freedom is excited (see Eq. (59)). In our new notation it can be written as 

X{t) = ljLl{t)'lpN/4 + lJi3{t)'tp3N/4: (64) 

where fxi{t) = fX3{t). Substituting (Xi{t) = fX3{t), fX2{t) — into Eqs. (63), we 
find that the first and third equations are both reduced to the equation of the 
harmonic oscillator 

ill + 2/ii = 0, (65) 

while both sides of the second equation turn out to be zero. Thus, 

(t) ^Acos{V2t + (t)o) (66) 

and we can linearize Eqs. (63) near the dynamical regime (66). Note that one 
can assume 0o = in Eq. (66), without loss of generality, since it is possible 
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to shift the time variable t in the dynamical equations of the bush B[a^] by 
the appropriate value. 



According to the general idea of linearization, let us suppose that 
l,,{t)^l,f\t)+5,it), 

li2{t) = 52{t), (67) 

Here 5i {i = 1, 2, 3) are small corrections to the original dynamical regime (64) 
^fii'\t), 0, /ii°''(t)| which is determined by the one-dimensional bush B[a^, ai]. 
Let us now substitute (67) into Eqs. (63) and neglect all nonlinear in Si terms. 
Then we have: 

Si + 2Si = -SaAcos{V^t)S2, 

S^ + 4^2 = -8aA cos(A/2i) (Si - S3) , (68) 
4 + 2^3 = 8aAcos{V2t)S2. 



Subtracting the third equation from the first equation and introducing new 
variables 71 (t) = Si{t) — 53(t), 72 (^) = '^2(^), we obtain finally 

71 + 271 = -16aAcos{V2t)-f2, 

(69) 

72 + 472 = — 8Q;Acos(\/2t)7i. 

Namely for these equations, we must find analytically the nontrivial solution 
(71(^)7^0,72(^)7^0) 

Note that the same system can be obtained from the linearized equations for 

investigation of the stability of the bush B[a^, a^i] which is presented in Table 8 
for the FPU-q; chain. Indeed, let us assume j = in the last system from this 
table, keeping in mind the relation k = j^. Then we have 

i>o = 0, 

l'N/4 + ^ll/A^N/A = -^^/^(^)^Af/2(t), 



(70) 

l'N/2 + UjI[I2^N/2 = -^^//(i) [VN/A{t) + i^3JV/4(^)] , 



where cVj = 4sin^(^), while /i(t) = Acos{^/2t) is the solution to the equation 
/i + 2/x = 0. Summing the second and forth equations of the system (70) leads 
us again to Eqs. (69), where 71 (t) = VN/iit) + vsN/iit), 72 (t) = i'N/2{t). 
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Thus, the analysis of the stabihty loss by the bush ai] is reduced to the 
problem of parametric excitation of the nontrivial solution to Eqs. (69). On the 
other hand, it is well known that the onset of unstable motion for the Mathieu 
equation can be investigated with the aid of different asymptotic methods of 
nonlinear mechanics (averaging method, multiscale method, etc.). In [20], we 
use one of such a type methods, namely, the normalization procedure based 
on the Poincare-Dulac theorem, for studying the parametric excitation in the 
system (69). Using this procedure, we remove, step by step, the so-called 
nonresonance terms of a given order to higher orders in the small parameter 
(in our case, this is the root mode amplitude A) with the aid of the appropriate 
nonlinear transformations of the dynamical variables ji{t). In such a way, the 
normalized system of differential equations with respect to new variables 7j(t) 
is obtained and, as a rule, this system turns out to be simpler than the original 
system. 

In [20], it was shown that the normalized equations up to 0{A^), obtained 
from Eqs. (69), can be solved exactly. The corresponding solution for 7i(t) 
contains the growing exponent exp(\/2a^A'^t) which increases from zero for 
the arbitrary small value of the root mode amplitude A. For the old variable 
7i(t) this exponent, being decomposed into the power series with respect to 
A, produces the unremovable secular term of the form 

a'^A'^\/2tsm{V2t + S). 

From the above said, it is clear that the parametric excitation in the sys- 
tem (69) takes place for the arbitrary small value of the root mode amplitude 
A and, therefore, the bush B[d'^, ai] possesses the stability domain of zero size, 
as it has already been concluded directly from Fig. 9. 

5.4 Diagrams of stability for one- dimensional bushes in the FPU-j3 chain 
54.1 BushB[a^,i] 

The stabihty diagram corresponding to the bush B[a^, i] in the FPU-/3 chain is 
presented in Fig. 10. From this diagram, one can see that for any finite N and 
for sufficiently small root mode amplitude A, there exists, in our mechanical 
systems, only one mode i'N/2{t), which is the root (and single!) mode of the 
bush B[a?, i]. On the other hand, if we fix arbitrary small value of A and then 
begin to increase the number of particles N, such can be always found for 

The full approximate analytical solution to Eqs. (69) represents the power series 
for 7i(t) and 72 (i) with respect to the small parameter A. Since this solution is 
rather cumbersome, we do not present it here. 
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which some other modes, close to the mode i'N/2{'t), turn out to be excited 
(see the center of the diagram 10). In other words, the threshold of the loss of 
stability of the considered bush decreases with increasing A^. Thus, the bush 
B[a^,i] becomes unstable in the continuum limit — > oo. 

Another interesting fact can be revealed with the aid of the diagram 10. 
Namely, there are certain intervals of mode numbers on r.h.s. and l.h.s. of 
this diagram which cannot be excited because of the interactions with the 
mode i'N/2{t). Note that the both above mentioned facts were found and 
discussed, in rather different terms, in the paper [14] . 



5.4.2 BushB[a^,i] 

The form of the bush B[a^, i] in the configuration and modal spaces was given 
in Eq. (57) The corresponding stability diagram is shown in Fig. 11. From this 
diagram, we find that, for small amplitudes, the bush B[a^,i] is not stable in 
the continuum limit N ^ 00, but can be stable for finite N. For example, 
for A^ = 12 the stability takes place up to the value Ai of the root mode^^ 
amplitude A. Note that three horizontal dashed lines are depicted in Fig. 11: 
the lower line is described by the equation A = Ai, while the middle and 
upper lines correspond to the equations A — A^ and A — A^, respectively. 

For Ai < A < A2 all the modes, with the only exception of the modes with 
j = ^ = ^,j = ^ = Y j ~ ~ excited. Therefore, for 

A^ = 12, the bush with trivial symmetry is excited. Indeed, the excitation of 
the mode with j = ^ can lead, in principle, to excitation of all the modes 
withj = ^(A; = l,2,---,ll). 

It is very interesting that there exists also the upper value (A3) for the loss of 
stability of the bush B[a^,i]. Indeed, for A > A^ (see Fig. 11) this bush turns 
out to be stable for all values of the number (A^) of particles in the chain. 
Thus, the bush B[a^,i] becomes stable even in the continuum limit, for such 
values of the root mode amplitude A. 



The authors of this paper used the direct integration of the dynamical equations 
of the FPU-/3 chain in the modal space, while we use the Floquet method. 

This mode is a certain superposition, eiipj^f/^ + S2tp2N/3j of the conjugate modes 
V'iv/s and ip'2N/3- 

Note that only the results of a linear stability analysis are presented in our stabil- 
ity diagrams: wc find those modes, which arc excited because of the direct parametric 
excitation with the root mode of the considered bush. On the other hand, if some 
new modes appear in the mechanical system, we must also take into account their 
interactions with other modes (it is the problem of nonlinear stability analysis). 
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5.4-3 Bush B[a'^,cd] 



The description of the bush B[a^,ai] in the configuration and modal spaces 
has already been given by Eq. (59): 



The stability diagram for this bush for the FPU-/? chain is represented in 
Fig. 12. From this diagram, one can see that, for any finite N, the bush 
B[a^,ai] is stable for sufficiently small root mode am})litudcs A. For example, 
for = 12, this bush loses its stability only for A > Ai (see Fig. 12). For 
very large A^, the bush B[a'^, ad] practically also remains stable for sufficiently 
small A because the intervals of the numbers of the excited modes turn out 
to be very narrow (these intervals are situated near j — ^ and j — ^). In 
the continuum limit N — > oo, rigorously speaking, the bush B[a^, di] appears 
to be unstable, but its instability is weak. 

We have just discussed the stability of all one-dimensional bushes for the 
FPU-/3 chain which are induced by the dihedral symmetry group. Now let 
us consider the stability of three "additional" one-dimensional bushes whose 
existence is brought about by the evenness of the potential of the FPU- (3 
chain. 

5.4-4 Bush B[d^,iu\ 

The description of this bush in configuration and modal spaces looks as follows: 



The stability diagram of the bush B[a^,m] is given in Fig. 13. It seems sim- 
ilar to that of the above discussed bush B[a^,i] (see Fig. 11). Therefore, the 
discussion of the stability properties of the bush B[a^,m] can be practically 
reduced to that of the bush B[a^, i]. 

5.4-5 Bush B[d'^,iu] 

This bush can be described by the equation 

X{t) — {x, —X, —X, x\x, —X, —X, x\ . . .} — i^(t)t/'Ar/4. (72) 

The corresponding stability diagram is represented in Fig. 14 and looks as 
"rabbit ears". For root mode amplitudes A > Ai (see Fig. 14) this bush is 
stable for arbitrary N. It is interesting that for the case = 12 there are 



X{t) — {0, X, 0, —x\0, X, 0, —X 



. . . } = e3u{t) [lpN/4 + 'il^3N/4] - 



X{t) = {x, —2x, x\x 



, -2x, X\...} = v{t) [£2'0iV/3 - £l'02iV/3 ■ 



(71) 
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no excited modes, certainly, besides the mode tf^N/A (or the mode il^sN/i for 
the equivalent bush B[a'^,a^m] which looks in the modal space as X{t) = 
z/(t)'037V/4)- This is in agreement with Table 10 from which we find that this 
bush is stable, at least, up to the root mode amplitude A < 20. 



5.4-6 Bush B[a^,ai,a^u] 
This bush can be described as follows: 
X{t) = {0, X, X, 0, —X, —x\0, X, X, 0, —X, —x\ . 



..} = u{t) [eiV'jv/e + £2-05^/6 

(73) 



One of the distinguishing features of the stability diagram for this bush (see 
Fig. 15) is the presence of three white regions in the "black sea" of instability. 
They correspond to the second stabihty zones for some modes whose num- 
bers fall into these white regions. Another interesting feature of the stabihty 
diagram for the considered bush is the very narrow unstable intervals in the 
"white sea" of stability (see the black lines which resemble parabolas). We 
did not examine these unstable parabolas in detail, because they seem to be 
not very essential for consideration of the stability of the bush B[a^, ai, ct'u]. 
Indeed, they represent too narrow regions of unstable movements. 



5.5 About numerical values of root mode amplitudes for the FPU chains 



The FPU models of a and (3 type, considered as exact mechanical systems, do 
not contain the interparticle distance a. As a consequence, one cannot estimate 
the reasonable values of the root mode amplitudes for the above discussed 
bushes. This is an important problem because atomic displacements in real 
crystals cannot, as a rule, exceed 10 — 15% of the equilibrium interparticle 
distance. 

On the other hand, we can consider FPU system as a certain approximation 

to more realistic models of one-dimensional crystals. Let us consider a nonlin- 
ear monoatomic chain with the Lennard- Jones potential U{r) describing the 
interaction between the nearest particles: 




Here, r is the interparticle distance, while e and a are certain constants. 

Let us expand the force /(r) acting on one particle, in Taylor series near the 
equilibrium interparticle distance a — a-\f2 up to the terms of the second 
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order: 



Then we can write the Newton equations of motion in the form 



(7^ ZCT"^ L 



(76) 



where m is the mass of one particle. 

Introducing new dimensionless time and space variables 



, 9\/4e , 21^ 

t = \ = Xn, (77) 

V ma^ a 

we obtain the following equations: 

x[ = (-2x; + x[^, + xl,) + [{x[^, - x[f - {x[ - xUn (78) 

These equations precisely correspond to the FPU-o; model considered in the 
present paper. 

On the other hand, the dimensionless variables x^, according to Eq. (77), 
are proportional to the ratio of the real atomic displacements Xi and the 
equilibrium interparticle distance a. As it has already been noted, for real 
crystal, Xi/a must not exceed 10 — 15%. The dimensionless displacements x'^ 
from the FPU-« model (78) can be considerably larger than this value, because 
of the coefficient 21 ~ 29.7 before Xi/a. For example, if Xi ~ 0.1a, then 

Let us now recall that, studying the bush stability, we use such normalization 
of the normal modes, that their amplitudes are of the same order as the atomic 
displacements corresponding to them. Thus, the root mode amplitudes A at 
about 3 seem to be rather reasonable. It is precisely because of this fact that 
we depict the bush stability diagram up to A of such an order of magnitude. 



6 Stability of two-dimensional bushes for the FPU chains 



In the previous section, we considered the stability of one-dimensional bushes 
with respect to increasing their root mode amplitudes. Such an increasing leads 
to strengthening of the parametric interactions of the root mode of a given 
bush with all the other modes which do not belong to this bush and, finally, to 
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a loss of its stability. Actually, the threshold of the loss of stability of any one- 
dimensional bush depends only on its total energy (the energy of the initial 
excitation), but not on the distribution of this energy between the kinetic 
and potential counterparts [20]. In two last columns of Table 10 we give the 
energy thresholds of the loss of stability for all one-dimensional bushes in the 
FPU chains. In contrast to one-dimensional bushes, the stability of the multi- 
dimensional bushes depends not only on the energy of the initial excitation, 
but also on the complete set of the initial conditions. Let us consider this 
problem in more detail. 

The four-dimensional phase space corresponds to the two-dimensional bushes 
considered below in the present paper. In Fig. 18-27, we depict only certain 
planar sections of the four- dimensional stability domains, aiming to show how 
nontrivial the boundary of these domains can be. 

The complete set of initial conditions for the dynamical equations of a given 
two-dimensional bush determines the concrete values of i/i(0), i'i(O), 1^2(0), ^'2(0), 
where and z>j(t) (i = 1,2) are the two modes of the considered bush 
and their velocities, respectively. We depict the bush stability diagrams (see 
Fig. 18-33) in certain planes which are chosen by fixing two of the four above- 
mentioned initial values. Most frequently, we specify 2^1(0) = 0, z>2(0) = and 
change i^i(O), ^'2(0) in some intervals near their zero values. Therefore, the ini- 
tial kinetic energy assumed to be equal to zero and the dependence of bush 
stability on the initial mode amplitudes i^i(O) and 1^2(0) is studied. 

The diagrams discussed in the previous section allow us to analyze stability of 
one-dimensional bushes for the FPU chains with arbitrary number (A^) of par- 
ticles. Unfortunately, we cannot obtain similar diagrams for two-dimensional 
bushes. Stabihty diagrams considered in this section depend on the number of 
particles and, therefore, they must be constructed individually for each spe- 
cific N. Because of this reason, we consider stability of 2D bushes only for 
fixed (for the most diagrams A^ = 12). As a consequence, we do not discuss 
the stability properties of the 2D bushes with translation symmetry [a^], [a^] 
and [a^°] (the complete hst of 2D bushes for the FPU chains can be found in 
Appendix) . 

Firstly, let us consider the stability of the two-dimensional bush B[a'', i] in the 
FPU-CK chain with A" = 12 particles. The stability domain for this bush, in 
the plane of the initial values of its both modes [z^i(O), ^'2(0)], is presented in 
black color in the center of Fig. 16. This domain, resembling a beetle, is 

In stability diagrams for one-dimensional bushes, we use white color for stabil- 
ity regions and black color for instability regions. In contrast, for two-dimensional 
bushes, white color corresponds to instability domains, while black color is used for 
stability domains. It seems to us that this manner of depicting stability diagrams 
turns out to be more expressive. 
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depicted on the background of the energy level lines. The bush B[a^,i] loses 
its stability when we cross the boundary of the black region in any direction. 
Moreover, it can be shown that the considered bush transforms into different 
bushes of larger dimensions at different points of the boundary. Which bushes 
appear near the various parts of the boundary of the stability domain of the 
bush B[a^,i] will be discussed elsewhere. Now we merely want to emphasize 
that the loss of stability of this bush is not determined only by the energy 
of the initial excitation (the boundary of stability domain does not coincide 
with any energy level line!), but depends on the relation between the values 
of its both modes [i/i(0), 1/2(0)], as well. Thus, the stability domain of the 
bush B[a^,i] is strongly anisotropic: the hnear size of this domain in different 
directions, from zero point ^liff) — 0, 1^2(0) = 0, turns out to be different. 

The stability diagram of the bush B[a^,i] in the FPU-a chain, depicted in 
Fig. 16, was obtained by the direct numerical experiment for the case z>i(0) = 
0, z>2(0) = 0. Therefore, initial kinetic energy was equal to zero and the total 
excitation energy was purely potential. 

It is interesting to consider the case where the excitation energy is distributed 
at the initial instant, in some a manner, between its kinetic and potential 
counterparts. The appropriate stability diagram is depicted in Fig. 17. Ac- 
tually, this is simply another section of the same four-dimensional stabil- 
ity domain of the considered bush which is determined by the equations 
z>i(0) = 0.2,z>2(0) =0.1. 

Comparing Fig. 16 and Fig. 17, one can see the effect of the "wind blowing 
from the right to the left". This effect leads to a loss of the symmetry of the 
diagram in Fig. 16 with respect to the vertical axis passing through the origin. 

Thus, the threshold of the loss of stability of the bush B[a^,i] in the FPU-a 
chain depend essentially not only on the distribution of the initial energy be- 
tween its two modes, but also between the kinetic and potential counterparts. 
Let us repeat, once more, that instead of dealing with the four-dimensional 
stability diagram, we investigate different sections of this domain. 

The stabihty diagrams for two other two-dimensional bushes, B[a^] and B[a^, ai], 
in the FPU-a chain are depicted in Fig. 18 and Fig. 20, respectively. These 
diagrams, as well as the diagrams in Fig. 21-27 for the FPU-/9 chain, corre- 
spond to the initial conditions z^i(O) = 0, z>2(0) = 0, i.e. the kinetic energy of 
the initial excitation is supposed to be equal to zero. 

Let us emphasize that the stability diagrams for the same bushes in the FPU-a 
and FPU-/3 chains look very different. As an example, one can compare Fig. 19 
and Fig. 22 for the bush B[a'', i] in FPU-o; and FPU-/3 chains, respectively. It is 
interesting to note that the last diagram (Fig. 22) demonstrates the presence 
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of several regions of instability (white color) in the stable sea (black color). 

In Figs. 24-27, we represent the stability diagrams corresponding to four ad- 
ditional two-dimensional bushes for the FPU-/3 chain (as it has already been 
said, they are induced by the evenness of the potential of this model). 

The most nontrivial picture corresponds to the stability diagram of the bush 
B[a^, a'^u]. It resembles a chivalrous cross with two crossed swords (along the 
diagonals of the diagram). Note that these swords describe, actually, the sta- 
bility of two equivalent one- dimensional bushes h'(t)[if)3+if)g] and i/(i)[i/?3— i/^g], 
while the two-dimensional bush B[a^, a^u] is i'i{t)i/}3 + i'2{t)i{)g. 

For many bushes, the stability domains depend essentially on the number N of 
particles in the chain. The area of the corresponding regions of stability, in our 

diagrams, tends often to decrease with the increasing of N. We demonstrate 
such a behavior in the cases N = 8,12,16,20,24,48 for the bush B[a'^,i] in 
the FPU-a chain (see Figs. 28-33) . Note that the stability domain of this 
bush shrinks only along the horizontal axis i^i(O) when N ^ oo. This is the 
consequence of the fact that when ui{t) — > the bush B[a^, i] — i^i(t)i/'3Ar/4 — 
J^2{t)il^N/2 (or its equivalent form B[a'^, ci^i] = z/i(t)'0jv/4 + i^2{t)'^N/2) tends to 
the one-dimensional bush B[a^,i] = i'2{t)'ipN/2- On the other hand, the latter 
bush possesses, for N oo, a. finite one-dimensional interval of stability with 
respect to increasing of the initial value ^'2(0) of its single mode 1^2 (^)- 



7 Conclusion 

The bushes of normal modes as exact nonlinear excitations in the physical 
systems with discrete (point or space) symmetry groups were introduced and 
discussed in our previous papers [2,3,4,6,9,10,11]. In general case, bushes can 
be rather complex, for example, the bushes in the Cgo fuUerene structure [11]. 

In [1], we considered bushes in monoatomic nonlinear chains which represent 
the simplest systems with translational symmetry. Mainly, we studied there 
the bushes induced by the group T of pure translations and discussed only 
in part those induced by the dihedral symmetry group D. The stability of 
a few simplest bushes was also studied in [1]. Nevertheless, there are some 
more onc-dimensional and two-dimensional bushes in the FPU-chains, in par- 
ticular, those obtained by Bob Rink in [15], whose stability still has not been 
investigated. On the other hand, the problem of the bush stabihty is one of 
the most important problems in the general theory of the bushes of normal 

We do not discuss the continuum limit — > 00 for two-dimensional bushes in 
the present paper. 
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modes. Because of this, in the present paper, we discuss in detail the stabihty 
of all one- dimensional and some two-dimensional bushes in both FPU-a and 
FPU-/5 chains. 

Let us summarize the results obtained in the previous sections of this paper. 

1. A simple crystallographic method for finding bushes in nonlinear chains is 
developed. 

2. The stability of all one-dimensional bushes in the FPU-« (three bushes) 
and in the FPU-/3 (six bushes, including those obtained in [15]) are inves- 
tigated. The stability diagrams which can be used for studying stability of 
one-dimensional bushes for the FPU chains with any finite and infinite num- 
ber of particles (continuum limit N — > oo) are presented. They were obtained 
by Floquet method. 

3. The stability diagrams of different type, which represent certain sections 
of four-dimensional stability domains, corresponding to some two-dimensional 
bushes in the FPU-a (three bushes) and in the FPU-/? (seven bushes), are 
found. These diagrams demonstrate explicitly that the bush stability domains 
depend not only on the energy of the initial excitation, but on the complete 
set of the initial conditions in the bush phase space. 

The above mentioned diagrams for one-dimensional and two-dimensional bushes 
show how nontrivial the shape of the stability domains for the bushes of nor- 
mal modes can be. 

In this paper, we consider only the main mechanism for the loss of bush 
stability, wherein it is brought about by the parametric excitations of the bush 
modes with those modes which do not belong to the given bush. Nevertheless, 
one can consider some different reasons for a bush to lose its stability, such 
as a finite temperature of the thermal bath, the presence of impurities in the 
chain structure, etc. These problems, as well as the methods for excitation of 
the bushes of normal modes in the FPU-chains will be considered in a separate 
paper. 
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Appendix. Two-dimensional bushes of vibrational modes for the 
FPU chains 



Table 6 

Representation of 2D bushes of vibrational modes for the FPU chains in the con- 
figuration space (classification by the dihedral group D) 



Bush 


Displacement pattern 




\X1,X2,X2,\ 


B[a^,i] 




T-k r A 9-1 

B[a^, a i\ 


\xi, -Xi,X2, -X2I 


B[a^i] 


|a;i,a;2,0,-X2,-xi| 


B[a^ ,ai] 


|0,xi,a;2, -X2, -xi\ 


B[a^a2^] 


|xi,-xi,a;2,0, -a;2| 




|xi,0, -Xl,X2,-2;2| 




.Cl. .('2. -Xo. -.1-1,0 


B[a^,ai] 


|0,xi,X2,0,-a;2,-a;i| 


B[a^,ct'i] 


Xi,0, -Xi,X2,0, -X2 


B[a^aS^] 


|xi, X2, 0, — X2, — Xl, 0| 
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Table 7 

Additional 2D bushes for the FPU-/? chain 



Bush 


Displacement pattern 


B[a^, aiu] 


\Xi,X2,X^,X2\ 


B[a^, a^iu] 


\x2,Xi,X2,Xz\ 


B[a^a2u] 


|xi,X2,-.xi,-a;2| 


B[a^,iu] 


|xi,X2,X3,X2,Xi| 


B[a^,aiu\ 


\XI,X2,X^,X^,X2\ 


B[a^, a^m] 


\xi,Xi,X2,X^,X2\ 


B[a^,a^iu\ 


\xi,X2,Xx,X^,X2,\ 


B[a^,aHu] 


\xi,X2,X2,Xi,X2,\ 


B[a^,iu] 


\xi,X2,X^,X^,X2,Xi\ 


B[a^, a^iu] 


\xi,Xi,X2,X2„X2„X2\ 


B[a^, a'^iu] 


|a;i,X2,.T2,.Xi,X3,X3| 


B[a^, i, a^u] 


\xi,X2,Xi, -Xi, -X2, -Xi\ 


B[a^,a'^i,a^u\ 


|xi, -Xi,X2,-Xi,Xi,-X2\ 


B[a'^,aH,a^u] 


\xi,X2,-X2,-Xi,-X2,X2\ 


B[a^, i, a'^u] 


|xi, X2, X2, Xi, -Xi, -X2, -X2, -Xi| 


B[a^,a'^i,a'^u\ 


|xi, -Xi, X2, X2, -Xi, Xi, -X2, -X2I 


B[a^,a^i,a'^u] 


|X1, X2, -X2, -Xl, -Xl, -X2, X2, Xl\ 


B[a^,a^i,a\] 


Xl,Xl,X2, -X2, -Xl, -Xl, -X2,X2 


B[a^, ai, a'^u] 


|0,Xi,X2,Xl,0, -Xl, -X2, -Xi| 


B[a^,a^i,a'^u\ 


Xl,0, -Xl,X2, -Xl,0,Xl, -X2 


B[a^,a^i,aSi\ 


|xi, X2, 0, -X2, -Xl, -X2, 0, X2I 


B[a^, aJi, a^u\ 


|X1, X2, Xl, 0, -Xl, -X2, -Xl, 0| 


B[a}^ ,ai,a^u] 


|0, Xl, X2, X2, Xl, 0, -Xl, -X2, -X2, -Xi| 


B[a}-^,a^i,a^u] 


|X1, 0, -Xl, X2, X2, -Xl, 0, Xl, -X2, -X2I 


B[a^^ ,a^i,a^u\ 


|xi,X2,0, -X2, -Xl, -Xl, -X2,0,X2,Xi| 


B[a}^ ,0^1,0.^ u] 


|xi,Xi,X2,0, -X2, -Xl, -Xl, -X2,0,X2| 


B[a}^,o?i,a^u] 


|X1, X2, X2, Xl, 0, -Xl, -X2, -X2, -Xl, 0| 
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Table 8 

Linearized dynamical equations for studying stability of ID bushes for the FPU-a 
chain. Here z>(t) is the vibrational mode representing a given bush, while k = ^ 
determines the number j of "sleeping" modes 



B[a'^,i]: i) + iiy = 



Uk + 4j^fc sin"^ k = — ^vu. 



. sin 2k 



viL-k + 4z^2t-fe cos^ k = —^uuk sm2k 



i>fc + Auk sin^ k 



— "VS t 



(1 - 2 cos {2k + f )) + (1 - 2 sin (2A; + f )) vi^^^ 

^l+fc + 4i/|+fc sin2 {k + f ) = (1-2 cos {2k + f )) i/fc + (1 + 2 cos 2k) 

i^2^+k + ^^if+k cos^ + f ) = [(1 - 2 sin {2k + f )) z/^ + (1 + 2cos2fc) 



i>k + sin2 k = [(1 - V2 - 

^>.+, + 4^.+,sin2 {k + l) = -^-^v 
i>|+, + 4z..+fcCos2A; = -^j/[(U 
f>^+fc + 4^^^+, cos^ (fe + f ) = 



cos (2A; + f )) V2L+k + (1 - V2sin {2k + f )) 1/^+^^ 
(1 - V2COS {2k + 5 )) z/fe + (1 + V2sin {2k + f )) 
(1 + x/2sin (2A; + f )) v^+k + (l + a/2 cos (2A; + f )) 1^3. 
(1 - V2sin {2k + I)) i/fc + (1 + V2 



cos {2k + l))v^ 
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Table 9 

Linearized dynamical equations for studying stability of ID bushes for the FPU-/3 

chain 



N ^ 



ui, + 1 ( 1 + u,, siu^ k = 



B[a^,a'^i]: i) 



i>k + 4:h'k sin^ k = — ^u^ sink 2vk sinfe — J^|+fc sin {k + ^) + fi^^j^j^ sin [k + ^) 

+4i/|+fcsin2 I) = -^^u'2^\ri{k + l) -z/^ sin A: + 2i/|+fc sin (/c + f ) - z/^ 
^^+k + + T") = -i^i^^sin {k + ^) t'fcsinfc - sin + f ) + 2i/s 



viL+k sin (A; + f) 



+ '^2._^;.sin {k+^) 
Uj. sin A; + 2i/| 



i+fcSin(/c + f) 



+ f)+2i^^^,sin (fc+f)" 



[vk + 4.{l + ^4v'')ukSiT?k = {) 



i^fc + 4j^fc sin^ k = — ^^v'^ sin A; 



2z/fc sin A: + J^ii+fc sin (A; + ^) — J^a^r _|_^ sin (A; + ^) 



27r> 
3 ; 



t'l+fc + 4z/2:+fe sin^ + f ) = ^""^ + f ) sin A; + 2z/2r_|_j^. sin (fc + |) + z^a^r , j^. sin (A; + 

3 3 ^ _ 3 3 

i^ai^j^. + ^y^jL+k ^i"^^ + T") = ^™ + T") ^^"^ ^ + ^^"^ + f ) + 2^22L+fe sin (A; + ^) 



B[a4,m]: i> + 2zy = 

i^fc + 4 (l + f z^fcsin^ A; = -^u'^u^_ksm2k 
viL-k + 4 (l + ^1/2) viL-kCos^k = -^fi^^L'ksm2k 



B[a^,a^i,a^u\: i> + i/ 

i>fc + 4^'jfc sin^ k = — ^^'^ sin A; 2^^ sin A; + sin (A; + |) — z^27r_,.j^. sin (A; + ^) 
i>|+fc + 4z/|+feSin^ (^ + f) = -^z^^sin(A; + |) sin A; + 2z/|+;t sin (A; + |) + sin (A; + ^) 

+ '^^■i2L+k sin^ + = -iv sin (A; + [-z/fe sin A; + un+k sin + f ) + '^^22L+k + ¥) 



Z^27T 
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Table 10 

Thresholds of bush stability for the FPU-a and FPU-/? chains with N = 12. Here 
we give the maximal value of the displacement x, as well as the maximal energy for 
which the corresponding bushes lose their stability 



Bush 


Displacement pattern 


Threshold of the bush stability 






in displacement 


in energy (per atom) 






FPU-a 


FPU-/? 


FPU-a 


FPU-/? 


B[a2,i] 




0.302 


0.112 


0.183 


0.025 


B[a^^] 


\x, 0, —x\ 


0.203 


0.268 


0.047 


0.085 


B[a^, ai] 


|0, a;,0, —x\ 





1.161 





0.675 


B[a^,i,a'^u] 


1 X ^ X ^ t)C ij X' 1 




> 20 




> 10000 


B[a^, ai, a^u] 


1 J ^ X' J X 5 ^ X 1 




0.488 




0.079 


B[a^,iu] 


1 X ^ ^X ^ X 1 




0.157 




0.074 
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Fig. 7. Regions of stability (white color) of different modes of the FPU-a chain, 
interacting parametrically with the one-dimensional bush B[a^,i]. 









Fig. 9. Regions of stability (white color) of different modes of the FPU-a chain, 
interacting parametrically with the one-dimensional bush B[a^,ai]. 
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Fig. 10. Regions of stability (white color) of different modes of the FPU-/? chain, 
interacting parametrically with the one-dimensional bush B[a^,i]. 
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Fig. 12. Regions of stability (white color) of different modes of the FPU-/? chain, 
interacting parametrically with the one-dimensional bush B[a^,ai]. 
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Fig. 13. Regions of stability (white color) of different modes of the FPU-/? chain, 
interacting parametrically with the one-dimensional bush B[a^,m]. 
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Fig. 14. Regions of stability (white color) of different modes of the FPU-/? chain, 
interacting parametrically with the one-dimensional bush B[a^,m]. 



21 



T 



T 



T 








TT 



Fig. 15. Regions of stability (white color) of different modes of the FPU-/? chain, 
interacting parametrically with the one-dimensional bush B[a^, az, a^u]. 
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Fig. 16. Stability diagram of the bush B[a*^, i] in the FPU-a chain with A'^ = 12 for 
the case z>i(0) = 0, z>2(0) = 0. 
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Fig. 17. Stabihty diagram of the bush B[a^, i] in the FPU-a chain with A'' = 12 for 
the case z>i(0) = 0.2, z>2(0) = 0.1. 
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Fig. 19. Stability diagram for the bush Fig. 22. Stabihty diagram for the bush 
B[o^, i] in the FPU-a chain. B[a^, i] in the FPU-/? chain. 




-0.5-0.3 "^^-^8 

Fig. 20. Stabihty diagram for the bush Fig. 23. Stabihty diagram for the bush 
B[a^, ai] in the FPU-a chain. B[a^, ai] in the FPU-/? chain. 
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Fig. 24. Stability diagram for the bush 
B[o^,a^'u] in the FPU-/? chain. 




Fig. 25. Stabihty diagram for the bush 
B[a^,am] in the FPU-/? chain. 
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Fig. 26. Stabihty diagram for the bush 
B[a^,iu] in the FPU-/? chain. 
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Fig. 27. Stabihty diagram for the bush 
B[a^,i,a^u\ in the FPU-^ chain. 
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Fig. 28. Bush B[a^,z] in FPU-o; with 
= 8. 








Fig. 29. Bush B[a'^,i] in FPU-a with 
= 12. 



Fig. 30. Bush B[a'^,i] in FPU-a with 
N = 16. 





Fig. 31. Bush B[a^,z] in FPU-a with 
N = 20. 




Fig. 32. Bush B[a^,i] in FPU-a with 
N = 24. 








Fig. 33. Bush B[a^^] in FPU-a with 
N = 48. 
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